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ABSTRACT 


Galerkin's  method  was  used  to  calculate  the  natural 
frequencies  of  vibration  of  tapered  cantilever  beams.  It 
was  found  that  the  values  obtained  from  Galerkin's  method 
converged  very  quickly,  and  for  small  tapers,  reliable  values 
for  the  fundamental  and  second  natural  frequencies  were  ob¬ 
tained  when  two  and  three  terms  were  used  in  the  trial 
solution . 

The  trial  solution  which  was  used  in  this  thesis  was 
found  to  be  unsuitable  for  calculating  natural  frequencies 
of  modes  higher  than  the  third  mode,  and  for  calculating 
the  fundamental  frequency  of  vibration  when  the  effects  of 
shearing  force  and  rotatory  inertia  were  considered. 

Tests  were  conducted  to  determine  the  natural  frequencies 
of  vibration,  and  it  was  found  that  the  experimental  fre¬ 
quencies  were  always  lower  than  the  theoretical  frequencies. 
The  largest  differences were  found  for  the  shorter  beams, 
and  the  beams  with  the  largest  thickness  taper. 

From  the  theoretical  values  which  were  obtained  for 
the  second  natural  frequency,  it  was  noted  that  for  small 
tapers,  the  frequency  of  vibration  decreased  almost  linearly, 
as  the  thickness  taper  was  increased.  This  fact  was  used  to 
propose  a  method  for  finding  the  second  natural  frequency  of 
vibration  of  a  beam,  by  using  a  linear  extrapolation. 
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NOTATION 


y  beam  deflection 

x  distance  along  the  beam  from  the  fixed  end 

t  length  of  the  beam 


a,  p  width  and  thickness  parameters 
E  Young ' s  modulus 

G  shear  modulus 

p  mass  density 

I  moment  of  inertia  of  the  cross  section 

A  area  of  the  cross  section 

I  r 

k  numerical  factor  which  depends  on  the  shape  of  the 

cross  section 

a)  frequency  of  vibration  ,  cps 

p  fundamental  frequency  of  vibration  of  a  uniform  beam,  cps 

y^  trial  solution  which  contains  i  terms 
til 

^  1  shape  function 


Subscript  o  denotes  values  at  the  clamped  end  of  the  beam. 
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CHAPTER  I 


INTRODUCTION 


1.1  Aim  of  Thesis 


This  thesis  is  concerned  with  determining  the  natural 
frequencies  of  vibration  of  straight,  tapered  cantilever 
beams  of  rectangular  cross  section.  Galerkin's  method  9,10 
will  be  used  to  determine  the  first  three  natural  frequencies 
of  vibration  for  linearly  tapered  cantilever  beams  when  the 
effects  of  shearing  force,  rotatory  inertia  and  damping  are 
neglected.  Galerkin's  method  will  also  be  used  to  find  the 
fundamental  frequency  of  vibration  of  a  beam  of  constant 
width  and  linearly  varying  thickness  when  the  effects  of 
shearing  force  and  rotatory  inertia  are  considered.  Finally 
this  method  will  be  applied  to  the  problem  of  finding  the 
fundamental  frequency  of  vibration  of  a  cantilever  beam  with 
nonlinear  tapers,  when  the  effects  of  shearing  force,  rota¬ 
tory  inertia  and  damping  are  neglected. 

For  the  special  cases  of  a  uniform  beam  and  linearly 
tapered  beams  which  are  tapered  to  zero  in  one  or  both  dir¬ 
ections  about  the  axis  of  symmetry,  closed  form  solutions 


have  been  found 


1,2,3 


The  frequencies  obtained  from 


*  Numbers  in  brackets  denote  references  contained  in  the 
Bibliography. 
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Galerkin's  method  will  be  compared  to  the  frequencies  ob¬ 
tained  from  the  closed  form  solutions  and  from  Martin's 


method 


8 


for  these  special  cases. 

Experiments  will  be  conducted  with  several  tapered 
cantilever  beams  to  study  their  behavior  when  they  are  sub¬ 
jected  to  a  cyclic  disturbance.  The  data  obtained  from 
these  experiments  will  be  used  to  plot  frequency  response 
curves  from  which  experimental  values  for  the  natural  fre¬ 
quencies  of  vibration  will  be  obtained. 


1.2  Historical  Introduction 


The  equation  which  governs  the  motion  of  a  vibrating 
beam  when  the  effects  of  shearing  force,  rotatory  inertia 
and  damping  are  neglected  is 


(1.2.1) 


When  considering  a  prismatic  beam, 
to 


El 


this  equation  simplifies 


> 


which  is  easily  solved  for  the  natural  frequencies  of  vibration 


when  the  boundary  conditions  are  applied 


For  the  case 
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of  a  cantilever  beam  the  boundary  conditions  are: 


y(o,t)  -  o 

y  (o,t)  =  o 

y8“(t,t)  «  0 

y  U,t)  =  0. 


In  187 7 j  Rayleigh  considered  the  effect  of  rotatory 
inertia  and  added  a  term  to  the  right  hand  side  of  equation 
(1.2.1)  to  account  for  the  added  inertia  of  the  beam,  and 
in  1921  Timoshenko  considered  the  effect  of  a  shearing  force 
on  the  motion  of  a  vibrating  beam.  When  these  two  effects 
are  considered^  the  equation  governing  the  motion  of  a 
prismatic  beam  becomes 
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One  of  the  earliest  methods  used  to  determine  the 
natural  frequency  of  a  tapered  cantilever  beam  was  applied 
by  Kirchoff  in  1879 ,  when  the  problem  of  a  vibrating  cone 
with  a  clamped  base  and  free  tip  was  considered,  Kirchoff 
found  an  infinite  series  solution  for  the  resulting  differ¬ 
ential  equation ,  which  gave  the  fundamental  frequency  of 
vibration  of  the  cone. 


In  1913  Ward 


extended  Kirchoff 1 s  method  to  include 


bars  of  rectangular  cross  section.  Ward  considered  the  case 
of  a  beam  with  an  area  of  zero  at  the  free  end,  and  used 
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the  free  end  of  the  beam  as  the  origin  of  the  co-ordinates . 
(Figure  1-1) 


y 

FIGo  1-1  WARD " S  BEAM 


Ward  let  the  width  vary  as  xn  and  the  thickness  vary  as  xras 
so  that  the  area  and  moment  of  inertia  of  the  cross  section 
could  be  expressed  as 
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-n+m 

~n+3m 
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where  and  are  the  area  and  moment  of  inertia  of  the 
cross  section  at  x  5  1,  By  separating  the  variables  in  the 
differential  equation 
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and  using  the  substitutions  (for  the  case  of  m 
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4x 
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Ward  found 
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(n+.l) 


=  0, 


(1.2.2) 


to  be  the  equation  of  motion  for  the  beam.  Equation  (1.2.2] 
is  a  Bessel  equation  of  order  (rrfl) ,  and  its  roots  give  the 
natural  frequencies  of  vibration. 


Using  essentially  the  same  approach,  Wrinch 
arrived  at  the  same  results  as  Ward  and  Kirchoff 


in  1922 


In  1962,  Housner  and  Knight ley  4  studied  the  vibra¬ 
tions  of  a  linearly  tapered  cantilever  beam,  by  using  a 
stepped  beam.  They  considered  the  beam  to  be  made  up  of  a 
finite  number  of  prismatic  sections,  and  wrote  the  following 
equations  for  the  shear,  bending  moment,  slope,  and  deflec¬ 
tion  for  each  section. 
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To  determine  the  natural  frequency  from  these  equations, 
Housner  and  Knight ley  started  with  an  approximate  mode  shape 
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and  used  successive  iterations  to  improve  the  mode  shape 
until  all  of  the  boundary  conditions  were  satisfied.  When 
the  boundary  conditions  were  satisfied,  the  corresponding 
value  of  cd  was  calculated. 


and  Lindberg 


modi* 


In  1963,  Leckie  and  Lindberg 
fied  the  method  of  using  lumped  parameters  to  determine  the 
frequencies  of  vibration  of  linearly  tapered  cantilever 
beams.  By  assuming  the  deflected  shape  of  each  element  to 
be  a  polynomial  of  the  form 

2  3 

y  =  a  +  bx  +  cx  +  dx  , 

Leckie  and  Lindberg  were  able  to  write  a  dynamic  stiffness 
matrix  which  related  the  shear  and  bending  moment  to  the 
slopes  and  deflections  of  the  ends  of  the  elements.  When 
the  boundary  conditions  were  satisfied,  the  parameters  a, 
b,  c,  and  d  could  be  determined  and  the  natural  frequencies 
calculated. 

An  attempt  to  apply  finite  difference  equations  to  the 
differential  equations  that  govern  the  motion  of  a  linearly 
tapered  cantilever  beam  when  the  effects  of  shearing  force 
and  rotatory  inertia  are  considered  was  made  by  Rissone  and 


Williams 


in  1965.  By  replacing  the  derivatives  in  the 


equations  of  motion 
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by  the  corresponding  finite  difference  expressions  and 
adding  the  equations  together,  Rissone^  and  Williams  were 
able  to  write  one  equation  in  the  two  variables  Y  and  tjj 
for  each  element  of  the  beam.  For  a  beam  made  up  of  n 
elements  this  procedure  gave  n  simultaneous  equations  which 
were  solved  for  the  natural  frequencies  of  vibration. 

Another  approach  to  the  problem  of  determining  the 
natural  frequencies  of  a  tapered  beam  is  to  treat  the  beam 


8 


considered  a  linearly 


as  a  single  unit.  In  1956,  Martin 
tapered  cantilever  beam  with  the  origin  of  the  co-ordinate 
system  at  the  clamped  end,  and  expressed  the  area  and 
moment  of  inertia  of  the  cross  section  as 


A  =  AQ  f (x,  a,  p) , 

i  -  i0  g(*s  P) , 

where  a  and  (3  are  the  width  and  thickness  parameters,  and 
Aq  and  are  the  area  and  moment  of  inertia  at  the  clamped 
end.  The  equation  of  motion  then  becomes 


where 
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Martin  assumed  that 


A  f(x,  a,  p)  y  , 


A  -  Aq  +  A-^a  +  A?p  +  A^ap  +  A^az  +  A^P^  +  ...  , 

2  2 
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where  A^  1  s  are  parameters  to  be  determined  and  1 s  are 
functions  of  x.  By  considering  beam's  with  small  tapers, 

(a  ^  0.5,  0.5)  terms  of  greater  than  second  order  in 

a  and  (3  become  small  and  were  neglected.  Substituting  the 
functions  for  A  and  y  into  the  differential  equation  for 
motion,  and  comparing  like  powers  of  a  and  (3  a  series  of 
differential  equations  was  formed,  which  could  be  solved 
for  the  values  of  A^ .  Using  these  values  of  A^  the  natural 
frequencies  of  vibration  were  calculated.  This  method  gave 
a  lower  bound  for  the  natural  frequencies  of  vibration. 


In  1965,  Rao 


found  that  by  applying  Galerkin ' s 


method  to  the  equation  of  motion  for  a  vibrating  beam, 
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acceptable  values  could  be  obtained  for  the  fundamental 
frequency  of  linearly  tapered  beams.  Rao ' s  values  were 
an  upper  bound  to  the  actual  frequencies,  and  were  obtained 
by  solving  only  one  equation.  The  details  of  this  method 
are  given  in  Chapter  II „ 


1 . 3  Summary 

The  problem  of  finding  the  natural  frequencies  of  vib~ 
ration  of  non-uniform  cantilever  beams  is  the  same  problem 
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as  determining  the  natural  frequencies  of  vibration  of 
turbine  blades.  This  problem  has  attracted  the  attention 
of  many  investigators 9  and  for  the  special  cases  of  a 
uniform  beam  and  linearly  tapered  beams  which  have  an 
area  of  zero  at  the  free  end,  closed  form  solutions  are 
available.  However,  in  most  applications  the  beams  are 
not  uniform,  nor  do  they  have  an  area  of  zero  at  the  free 
end.  When  this  is  the  case,  approximate  methods  have  to 
be  used  to  determine  the  natural  frequencies  of  vibration. 

When  the  problem  is  analyzed  using  a  finite  difference 
method,  of  the  type  used  by  Rissone  and  Williams,  or  a 
lumped  parameter  method,  of  the  type  used  by  Housner  and 
Knight ley,  and  Leckie  and  Lindberg,  each  element  of  the 
beam  must  be  considered  separately,  and  the  number  of  equa¬ 
tions  which  are  formed  is  equal  to  the  number  of  elements 
considered.  To  solve  these  equations,  it  is  usual  to  use 
a  digital  computer,  and  the  accuracy  of  the  results  depends 
on  the  number  of  elements  used.  For  tapered  beams  the  re¬ 
sults  converge  slowly,  and  as  a  result  a  .large  number  of 
elements  have  to  be  considered.  The  results  obtained  using 
these  methods  are  always  higher  than  the  experimental  values, 
and  the  values  calculated  using  the  closed  form  solutions. 

When  the  problem  is  considered  by  treating  the  beam 
as  a  continuous  unit,  as  done  by  Martin  and  Rao,  it  was  found 
that  the  number  of  equations  which  have  to  be  solved  in  order 
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to  determine  the  frequencies  is  much  less  than  the  number 
of  equations  which  have  to  be  solved  for  the  finite  dif¬ 
ference  and  lumped  parameter  methods,,  As  a  result,  a  desk 
calculator  can  be  used  instead  of  a  digital  computer .  In 
view  of  this  fact,  it  appears  that  these  methods  of  solution 
are  more  practical  than  the  finite  difference  and  lumped 
parameter  methods .  However,  on  surveying  the  literature 
available,  it  was  found  that  Galerkin 1 s  method  has  only  been 
used  to  find  the  fundamental  frequency  of  vibration  of  lin¬ 
early  tapered  cantilever  beams,  and  before  any  judgement 
can  be  made  on  the  actual  merit  of  Galerkin “ s  method  as  a 
means  of  determining  the  natural  frequencies  of  vibration, 
it  must  be  applied  to  a  wider  range  of  problems. 

In  Chapter  II,  Galerkin* s  method  is  used  to  determine 
the  first  three  natural  frequencies  of  vibration  of  linearly 
tapered  cantilever  beams.  It  is  also  applied  to  the  problem 
of  determining  the  fundamental  frequency  of  a  beam  which 
has  a  constant  width  and  a  linearly  varying  thickness,  when 
the  effects  of  shearing  force  and  rotatory  inertia  are  con¬ 
sidered.  Finally  Galerkin' s  method  is  used  to  determine  the 
fundamental  frequency  of  vibration  of  a  cantilever  beam  with 
nonlinear  tapers.  The  convergence  of  the  results  is  dis¬ 
cussed  and  the  results  are  compared  to  the  results  which  were 
obtained  from  the  closed  form  solution  and  Martin's  solution. 


CHAPTER  II 


THEORETICAL  ANALYSIS  AND  RESULTS 

2.1.  Plan  of  Chapter 

In  Section  2.2,  the  equations  of  motion  for  vibrating 
cantilever  beams  are  formulated  when  the  effects  of  taper, 
shearing  force  and  rotatory  inertia  are  considered.  The 
equations  for  a  uniform  beam  with  and  without  the  effects 
of  shearing  force  and  rotatory  inertia  are  presented  with 
their  closed  form  solutions  in  Subsection  2.2-1.  In  Sub¬ 
section  2.2-2  the  equations  of  motion  which  were  used  to 
determine  the  natural  frequencies  of  vibration  of  tapered 
beams  are  presented. 

In  Section  2.3  Galerkin ' s  method  is  used  to  calculate 
the  natural  frequencies  of  vibration  of  tapered  cantilever 
beams.  Subsection  2.3-1  contains  a  brief  discussion  of  the 
requirements  of  Galerkin ' s  method,  and  in  Subsection  2.3-2 
the  one,  two,  and  three  term  trial  solutions  which  were 
used  in  the  application  of  Galerkin 's  method  are  derived. 

The  three  trial  solutions  are  applied  to  the  equation  of 

ij 

motion  of  a  linearly  tapered  beam,  and  the  resulting  equation 
which  must  be  solved  to  determine  the  natural  frequencies  of 
vibration  are  given  in  Subsection  2.3-3.  The  one  term  trial 
solution  is  applied  to  the  equation  for  a  beam  of  constant 


‘ 

' 
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width  and  linearly  varying  thickness  when  the  effects  of 
shearing  force  and  rotatory  inertia  are  considered,  and 
to  the  equation  for  a  nonlinearly  tapered  beam  in  Subsections 
2.3-4  and  2.3-5  respectively. 

The  results  obtained  for  the  natural  frequencies  are 
presented  in  Section  2.4.  In  Subsection  2.4-1  the  conver¬ 
gence  of  Galerkin ' s  is  illustrated  by  comparing  the  frequ¬ 
encies  obtained  when  the  one,  two,  and  three  term  trial 
solutions  were  applied  to  the  equation  of  motion  for  a  lin¬ 
early  tapered  cantilever  beam  to  the  frequencies  obtained 
from  the  closed  form  solutions.  In  Subsection  2.4-2  the 
frequencies  obtained  from  Galerkin* s  method  for  the  first 
three  modes  of  linearly  tapered  cantilever  beams  are  pre¬ 
sented  and  discussed.  The  frequencies  obtained  when  shear¬ 
ing  force  and  rotatory  inertia  are  considered  are  presented 
in  Subsection  2.4-3,  and  the  fundamental  frequencies  of  beams 
with  nonlinear  tapers  are  given  in  Subsection  2.4-3. 

In  Section  2.5  the  theoretical  results  and  the  problems 
encountered  in  calculating  them  are  discussed. 

2.2  The  Equations  of  Motion 

In  general,  it  is  assumed  that  the  beam  material  is 
homogeneous  and  isotropic,  therefore  the  mass  density  p, 
Young's  modulus  E,  and  the  shear  modulus  G  are  constant. 
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The  clamped  end  of  the  beam  is  the  origin  of  the  co-ordinate 
system  and  the  positive  direction  is  taken  downward  (See 
Figure  2-1) .  At  the  origin,  the  area  and  moment  of  inertia 
of  the  cross  section  are  given  by  AQ  and  I„  respectively* 

The  effect  due  to  damping  is  neglected. 


FIG.  2-1  CO-ORDINATE  SYSTEM 


2,2-1  The  Uniform  Beam 

The  equation  of  motion  for  a  vibrating  uniform  beam 
can  be  shown  to  be 


(2,2a) 


when  the  effects  of  shearing  force  and  rotatory  inertia  are 
neglected.  When  the  clamped  end  of  the  beam  is  taken  at 
the  origin  of  the  co-ordinate  system,  the  boundary  conditions 
for  a  cantilever  beam  are 


Y(0,t) 


0 


Y  (0,t) 


0 


i  e 


Y  (*,t) 


0 


1  I  0 


Y  (*,t> 


0, 


*)  ’  y  v 
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where  primes  denote  differentiation  with  respect  to  x0 
By  assuming  Y  to  be  of  the  form 

Y(x,t)  ™  y(x)  sin  cat  , 

Equation  (2.2,1)  reduces  to 

EX  -  pA  a)  y  ,  (2.2.2) 

dx 

and  the  corresponding  boundary  conditions  are 


y ( o)  -  o 

y!(0)  -  0 

y  (i)  *  0 

0  0  8 

y  U)  =  0  . 


(2.2.3) 


The  general  solution  of  Equation  (2.2.2)  is 

y(x)  ~  cosh  (kx)  +  C.7  sinh  (kx)  +  CU  cos  (kx.) 

+  sin  (kx) , 

where  k4  35 

El 

When  the  boundary  conditions  (2.2,3)  are  satisfied  the  char¬ 
acteristic  frequency  equation  becomes 


cos  (k t)  cosh  (kt) 


-1 


(2.2.4) 


and  the  solution  of  the  equation  of  motion  becomes 


Y(x,t)  -  C0  (cosh  (kx)  -  COS  (kx))  -  T  ifnlfk 


(sinh  (kx)  -  sin  (kx)) 


sin  (out)  . 
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The  first  three  roots  of  Equation  (2.2.4)  are 


V 


V 


1.875 


4.694 


n. 


n. 


7.855  -  n. 


and  the  corresponding  natural  frequencies  of  vibration  ares 


l 


El 

pA 


(2.2.5) 


Shearing  force  and  rotatory  inertia  are  two  effects 
which  reduce  the  natural  frequency  of  vibration.  When  shear¬ 
ing  force  is  considered  in  the  calculation  of  beam  deflections, 
the  deflections  increase,  which  means  that  the  stiffness  of 
the  beam  is  reduced  and  there  is  a  corresponding  reduction 
in  the  natural  frequency.  Rotatory  inertia  is  present  be¬ 
cause,  as  the  elements  of  the  beam  translate  they  also  ro¬ 
tate,  and  in  the  case  of  a  vibrating  system  the  rotation  of 
the  beam  adds  inertia  to  the  system.  This  additional  inertia 
has  the  effect  of  increasing  the  mass  of  the  system  and  results 
in  a  further  reduction  of  the  natural  frequency  of  vibration. 

When  considering  the  combined  effects  of  shearing  force 
and  rotatory  inertia,  the  slope  of  the  beam  is  assumed  to 


consist  of  two  components 


1 


which  is  the  slope  when 


shearing  force  is  neglected,  and  £,  which  is  the  angle  of 


shear . 


. 
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dy  _ 
dx 

» 

The  inertial  moment  caused 
the  beam  is 


f  +  £  . 

by  the  angular  acceleration  of 


dx  . 


When  shearing  force  is  considered,  the  shear  on  the  face  of 
any  element  becomes 

v  «  k ' £  AG 

8  » 

where  k  is  a  numerical  factor  which  depends  on  the  shape 
of  the  cross  section . 

Considering  these  new  terms,  the  equation  for  the  ro¬ 
tation  of  the  beam  becomes 


AG 


p1 


b2* 

at2 


(2.2.6) 


and  the  equation  for  the  translation  of  the  beam  becomes 


(2.2,7) 


For  the  case  of  a  uniform  beam  ip  can  easily  be  eliminated 
from  Equations  (2.2.6)  and  (2.2.7),  leaving  the  following 
equation  in  y. 


El 


+ 


0 


(2.2.8) 


Solving  Equation  (2.2.8)  and  applying  the  boundary  conditions 
for  a  cantilever  beam  gives  the  following  frequency  equation. 


.. 


■ 
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(A4  +  B4)  +  2A2B2  cos  (At)  cosh  (Bt)+  (A3B  -  AB3) 

sin  (At)  sirih  (Bt)  88  0 


(2.2,9) 


where 


A 


a  +  Va2“4B 


2  2 

E  •  » 

k  G 


2  4 

p  0) 
k  EG 


pAo)2 

El 


Another  case  when  f  can  be  eliminated  from  Equations  (2.2.6) 
and  (2.2.7)  is  for  a  beam  of  constant  width  and  linearly 
varying  thickness.  This  case  is  discussed  in  Subsection 
(2.2-2) . 


2„2“2  The  Tapered  Beam 

,SV  ft 

Consider  the  beam  shown  in  Figure  (2-2) ,  which  has  a 
rectangular  cross  section  and  length  t.  Let  it  be  clamped 
at  the  large  end  and  free  at  the  small  end. 


FIG.  2-2  TAPERED  BEAM 
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Defining  the  taper  parameters  as 


a  - 


bQ  bt 

b0 


and 


P  = 


tQ  ~  tt 

t0 


the  width  and  thickness  of  the  beam  at  any  point  a  distance 
x  from  the  origin  is 


b  =  bQ(  1  -  a(j) 


n 


t0  (  1  ^  ^ 


m 


Defining  In  and  as  the  moment  of  inertia  and  the  area 
of  the  cross  section  at  the  clamped  end,  the  moment  of  inertia 
and  area  at  any  point  in  the  beam  may  be  expressed  as 


I  =  I  (l~~aXn)  (l-pX™)  3  , 


( 2  0  2 


A  =  Aq  (l-aXn)  ( l-fiy?1)  , 


(2.2.11) 


where  X 


x 

l  * 


For  a  non  uniform  beam,  the  equation  of  motion  is 


E 


0  9  8  0  SO 

I  Y  +  21  Y 


+  IY 


where  primes  denote  differentiation  with  respect  to  x  and 
dots  denote  differentiation  with  respect  to  t.  Using  the 
above  expressions  for  I  and  A  and  assuming  that  Y(X,t)  =  y(x) 
sin  cjot ,  the  equation  of  motion  for  a  tapered  beam  becomes 


. 
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V 


-n  (n-1)  a  Xn~2  ( l-pX™-)  3  +  SmnapX11^11"2  ( l-pX™)  2 
-  3m(m-l)  (3X™  "2  ( l-aXn)  (l-px™)2  +  6m2£2X2m“2 
( l-aXn)  (l-px™) 


|  naXn  1(l-pXm)3 


+  3mpXm  1(l-aXn)  (l-px™) 


(l-px™)3  y 


IV 


0)2pA 


I  I  I 


+  (1-aX11) 


El 


0 


“  (l~aXn) (l-pxm)y  .  (2.2.12) 


For  a  linear  taper,  n  and  m  are  both  equal  to  one,  and 
Equation  (2.2.12)  simplifies  to 


( 1-ctX)  ( 1-(3X)  3  yIV  -  | 


a(l-pX) 3  +  3P(l-aX) (1-px) 2 


V 


6(32(l-aX)  (l-f3X)  +  6aP(l-px)2 


2  A 

a)  pA 


0 


El 


0 


(2.2.13! 


2.3 


and 


(1-aX) (1-PX)  y. 

For  the  limiting  cases  of  a  wedge  (a  =  0,  (3=1) 
a  cone  (a  =  (3  =  1)  3  ,  Equation  (2.2.13)  can  be  reduced  to 

a  Bessel  equation,  which  can  be  solved  for  the  natural  fre¬ 
quencies  of  vibration.  However,  there  is  no  general  analy¬ 
tical  solution  to  Equation  (2.2.13). 

When  the  effects  of  shearing  force  and  rotatory  inertia 
are  considered  in  the  problem  of  determining  the  natural 
frequencies  of  vibration  of  tapered  cantilever  beams,  it 
becomes  very  difficult  to  eliminate  either  Y  or  f  from  the 
equations 


:  (  ^X^vnrac'  +  o(X-a) o- 

;X«:  (X-  m)  !*£ 


.  £3 
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El 


'X 


pA 


+  k 


d  Y 
bx 


-  f 


AG  -  pi  = 

at2 


0 


a2y  .  •  a 
^  “  k  ^ 


^  -  *)  ™ 


=  0  - 


(2.2.6) 


(2.2.7) 


* 


One  case  for  which  f  has  been  eliminated  is  for  a  wedge 
The  procedure  used  to  eliminate  f  is  outlined  in  the  follow¬ 
ing  discussion. 

Since 

=  t  +  i 


and  Y(x:,t)  =  y(x)  sin  (art)  ,  f(xyt)  can  be  assumed  to  be 


f{xyt)  =  f{x)  sin  (cot)  . 

Substituting  these  values  for  Y(x,t)  and  fix,  t)  into  Equations 
2.2.6  and  2.2.7  produces  the  two  coupled  ordinary  differen¬ 
tial  equations 


_d_ 

dx 


El 


d  jr 

dx 


+  k 


dy 

dx 


AG 


2t  , 

po)  If 


(2.2.14) 


d 

dx 


k  AG  (1^  ■  *) 


2 

pco  Ay 


(2.2.15) 


For  a  wedge,  (a  =  0)  the  area  and  moment  of  inertia  of  the 
cross  section  is  given  by 

A  =  AQ(1  -  PX) , 

I  =  IQ(1  -  (3X)3  , 

where  X  =  y  . 


* 


The  term  "wedge"  is  used  here  to  describe  a  beam  of  con¬ 
stant  width  (a  =  0)  and  linearly  varying  thickness  (p  =  constant) 


' 
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Substituting  these  values  for  A  and  I  into  Equations  (2.2.14) 
and  (2.2.15)  give  the  following  equations  for  the  rotation 
and  translation  of  the  wedge 


JL 

dx 


E  ( 1- (3X)  3 


+  D(l-pX) 


(|Y  _  y)  m  _  po)2(l-px)3^,  (2.2.16) 


_d_ 

dx 


(1-PX) ( 


dv 

dx 


where 


k  GA 


D  = 


0 


0 


2 

-  (l-pX)y  , 

k  G 


To  simplify  these  equations,  let 


(2.2.17) 


(1  -  0X)  «  u  , 
then  d/dx  =  -\$/l)  (d/'du)  . 

Substitute  u  into  Equations  (2.2.16)  and  (2.2.17)  and  dif¬ 
ferentiate  (2.2.16)  with  respect  to  u  to  get  the  following 


r.  2 ,  ,  3/i 

3u  ip  +  u  ip  f 

2 

(cy  +  ip)  +  u(cy  +  ip  )  =  -  ±L-~-  uy  , 

ck  G 

where  c  -  (f3/£)  ,  and  the  primes  represent  differentiation 
with  respect  to  u.  By  employing  a  series  of  differentiations 
and  substitutions,  ip  can  be  eliminated  and  the  resulting 
equation  of  motion  becomes 

E  ( l“f3X)  2  (D  -  pcjcP  ( 1“(3X)  2  +  Ec2) 


1  d  y 

4  4 

cl  dX 


0  2  8  8  3 

6u  Tp  +  6u  7p  +  u  ip 


S  B 


+ 


Ppm 

« 

a 


uy 


•pm 


' 

. 
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+ 


+ 


„  3 

3 

l 

dX 

1 

.2 

d y 

2  2 

ct 

dX 

1R 

+  Ec2) 

k  G 


E(1-PX)(6D  -  4pa)2  ( 1~(3X)  3  +  6Ec2) 

(D  +  Ec2)(6Ec2  +  pa)2  ( 1-(3X)  2)  -  p2cd4(1-PX)4 

Epob2  ( l-f3X)  2  -  (l-px)4 

k  G 


1  dy 
‘  l  dX 


(D  +  Ec2)  3 pa)2  ( l-f3X)  +  —  ft- (1-pX) 


k  G 


2  4 


3Ep  od*  , ,  3  2  4,.  _ 3 

-  — -  (1-px)  -  p  a)  ( 1-pX) 

k  G 


+  y 


(D  +  Ec  ) 

0 

ck  G 


3Ec2 poo2  -  Dpa)2  +  2P2cd4  (1~PX) 


2  4 


3  6 

^-T-  (1-PX) 
ck  G 


=  0 


(2.2 


The  details  of  the  derivation  of  this  equation  are  given 
in  Appendix  I-A, 

Equation  (2.2.18)  is  used  in  Subsection  2.3-4  to  de¬ 
termine  the  effects  of  shearing  force  and  rotatory  inertia 
on  fundamental  frequency  of  a  wedge. 


2.3  Determination  of  the  Natural  Frequencies  Using 
Galerkin 1 s  Method 

2.3-1  Galerkin 1 s  Method 

This  method  was  proposed  by  Galerkin  in  1915,  as  a 
means  of  finding  an  approximate  solution  to  boundary  value 


.18) 
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The  solution  of  the  homogeneous  differential  equation 


+  an(x)z  =  0 


can  be  approximated  by  the  partial  sum 

n 


z  (x)  =  >  c .  <t> .  (x) 

n  1  l  ' 

i  =1 


(2.3.1) 


when  the  functions  $^(x)  are  selected  to  satisfy  the  bound¬ 
ary  conditions  of  the  problem.  The  approximate  solution 
z^(x) ,  however  does  not  exactly  satisfy  the  differential 
equation  L(z)  =  0,  and  when  z^(x)  is  used  in  the  differen¬ 
tial  equation,  an  error  is  present.  By  selecting  the  con¬ 
stants  c.  such  that 

l 


l 


I 


L(zn)  ^(x)  dx  =  0  1  =  1,  2,  ...  n 


o 


the  error  is  minimized,  and  the  expression  (2.3.1)  becomes 
an  acceptable  approximation  to  the  solution. 

2.3-2  Determination  of  the  Trial  Solutions 

In  order  to  apply  Galerkin ' s  method  it  becomes  nec¬ 
essary  to  find  a  set  of  functions  <t>^(x,t)  such  that  the 
boundary  conditions  are  satisfied.  When  the  functions 
<$>^(x,t)  are  found,  the  motion  of  the  beam  can  be  approximated 


X 

■■  x 

v 

' 
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by 


Yn(x,t) 


n 


a±  4>i(x}t)  . 


i  -1 

Since  Y(xat)  “  y(x)  sin  (cot)  9  <t»^(x,t)  can  be  approxi¬ 
mated  by  4>^(x)  sin  (a)t)  .  The  boundary  conditions  which 
0^(x)  must  satisfy  are 


<t>±  (0)  =  0 

<^'(0)  -  0 

*±"W  =  0 

8  0  8 

0±  (^)  =  0  » 

Let  y(x)  be  a  polynomial  of  the  form 

y(x)  "  cQ  +  c-^X  +  c2X2  +  c^X3  +  .  ..  +  c^X-5  . 


When  the  boundary  conditions  are  satisfied  with 
y^(X)  is  found  to  be 


and 


y (X)  *  a-,  ( 6X2  -  4X3  +  X4)  , 
<t>:i  (X)  -  (6X2  -  4X3  +  X4}  , 


j 


( 2  „  3  0  2) 


for  j  ”  5S  y2(X)  becomes 

y2(X)  -  a1(6X2  -  4X3  +  X4)  +  a2(20X2  -  10X3  +  XD)  , 
d2(X)  -  (6X2  -  4X3  +  X4)  , 

0  (X)  -  ( 20X2  -  10X3  4-  X5)  ,  (2.3.3) 

and  for  j  =  6,  y3 (X)  becomes 

y  (X)  -  a1(6X2  -  4X3  +  X4)  +  a2(20X2  -  10X3  +  X5) 
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^ x (X)  =  (6X2  - 
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4X  +  X  )  , 

4>2(X)  =  (20X2  - 

-  iox3  +  x3)  , 

<t>3(X)  =  ( 45X‘i'  - 

~  20X3  +  x6)  , 

(2.3.4) 

where  X  =  —  . 

This  process  can  be  carried  out  indefinitely,  and  each  new 
trial  solution  for  the  deflected  form  is  slightly  better 
than  the  previous  approximation. 


2.3-3  Natural  Frequencies  of  Linearly  Tapered  Beams 


The  governing  differential  equation  for  a  linearly 
tapered  cantilever  beam  when  the  effects  of  shearing  force 
and  rotatory  inertia  are  neglected  is 
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where  primes  denote  differentiation  with  respect  to  x, 
and  X  =  —  . 

Using  the  one  term  trial  solution 

yl(X)  =  ax(6X2  ”  4x3  +  X4)  ,  (2.3.2) 

in  the  application  of  Galerkin ' s  method  to  Equation  (2.3.5) 
gives  the  fundamental  frequency  of  vibration  when 
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4>  ]  dX  =  0. 


After  the  substitutions  are  made  for  y-,  and  the  derivative: 
of  this  integral  becomes 
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When  this  integration  is  carried  out,  the  fundamental  freq¬ 
uency  of  vibration  is  found  to  be 
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When  the  two  term  trial  solution, 
y  (X)  =  a,  (6X2--4X3+X4)  +  a0  ( 20X2-10X3+X5)  , 


(2.3.3) 
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is  used  in  the  application  of  Galerkin ' s  method  to  Equation 
(2.3.5),  two  equations  are  formed. 
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After  the  integrations  are  completed,  the  equations  become 
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In  order  to  get  a  non-trivial  solution  to  these  two  equations 


the  determinant  of  the  coefficients  of  an  and  a0  must  be 
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set  equal  to  zero.  Writing  Equations  (2.3.6)  and  (2.3.7)  as 
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where 
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When  this  determinant  is  expanded,  and  set  equal  to  zero, 
the  frequency  equation  becomes 
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Equation  (2.3.8)  can  readily  be  solved  for  two  values  of 
2 

k.  from  which  the  first  two  natural  frequencies  of  vibration 
are  determined. 

When  the  three  term  trial  solution, 


y3(X)  =  ax (6X2-4X3+X4)  +  a2 ( 20X2-10X3+X5)  +  a, (45X2-20X3+x6) , 

(2,3.4) 

is  used  in  the  application  of  Galerkin's  method  to  Equation 
(2.3.5),  three  equations  are  formed. 
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The  determinant  of  the  coefficients  of  a^,  a2  and  a^  must 
be  set  equal  to  zero  to  get  a  non-trivial  solution  to 
Equations  (2.3.9),  (2.3.10),  and  (2.3.11).  The  coefficients 
of  a-^  and  a2  in  Equations  (2.3.9)  and  (2.3.10)  are  identical 
to  the  coefficients  of  a^  and  a ^  in  Equations  (2.3.6)  and 
(2.3.7).  If  the  remaining  coefficients  are  defined  as 
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Equations  (2.3.9),  (2.3.10),  and  (2.3.11)  can  be  expressed 
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Observing  that  A13  =  A31,  and  A23  =  A32  the  determinant  of 
the  coefficients  of  a^,  a2,  and  a3  becomes 
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In  order  to  determine  the  natural  frequencies,  the 

determinant  of  the  coefficients  must  be  set  equal  to  zero. 
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i.e.  |  A  -  /c  B  |  =  0. 

When  the  determinant  is  expanded,  it  becomes  necessary  to 

find  the  differences  between  very  large  numbers.  These 

differences  are  so  small  that  the  round  off  errors  are 

large  enough  to  produce  erroneous  results.  To  avoid  this, 
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the  values  of  /c  are  found  in  the  following  manner. 
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where  V  is  a  column  vector  consisting  of  the  elements  a^ 
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which  can  be  reduced  to 

•1 


[if 


-1 


V  =  /c2V  , 


and  the  problem  of  finding  the  values  of  k  becomes  the 
problem  of  finding  the  characteristic  roots  of  the  matrix 


-1 


A 


B 


iT 


This  was  done,  without  expanding  the  matrix,  by  using  an 


iterative  procedure 


11 


2.3-4  Wedges  With  Shearing  Force  and  Rotatory  Inertia 


The  equation  of  motion  for  a  wedge  when  the  effects 
of  shearing  force  and  rotatory  inertia  are  considered  is 
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+  y 


)  (/3Ec2pcD2  -  DpcD2  +  2p2CD4(l-pX)  2N) 
L  ck  G  V 


0 


(2,2.18) 


When  this  equation  is  multiplied  by  c,  and  the  coefficients 
of  the  powers  of  cd  are  collected  it  becomes 

A^  (1-pX)4  y  cd6 


A2(1-PX)4  y"  +  A3(1-=PX)3  y'  +  A4(1-PX)2  y 


CD 


&  TV  T  lit  pii 

A5(l-px)  yxv  +  A6(1-PX)J  y  +  A?U-pX)x  y  +  Ag(l-px)  y 


+  S  y 


2  . 
a)  + 


A10(1-PX)2  yIV  +  A1X(1-PX)  y  +  A,,y 


III  II 

4-  a 

12- 


=  0, 


(2.3,13) 


where  primes  denote  differentiation  with  respect  to  X,  and 
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t3 


When  the  one  term  trial  solution  for  y(x) , 

y,  (x)  *  a.,  (6X2  -  4X3  +  X4)  # 


(2.3 


is  used  in  applying  Galerkin ' s  method  to  Equation  (2.3.13) 

2 

a  cubic  equation  in  a)  is  produced.  (See  Appendix  I-B) , 

2.3-5  Fundamental  Frequency  of  Nonlinearly  Tapered  Beams 

For  beams  with  nonlinear  tapers,  where  the  moment  of 
inertia  and  area  can  be  expressed  as 


I  =  IQ  ( l-aX11)  (l-px”1) 


(2.2.10) 
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A  =  Aq  ( l~aXn)  (l-px™) 
the  equation  of  motion  is 


(2.2.11) 
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03  pA 

=  t  (l-axn)  (l-pX^y  (2.2.12) 

EI0 


Using  the  one  term  trial  solution  (y^(x))  in  the  application 
of  Galerkin 8 s  method  to  this  equation  gives 
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which  can  be  solved  for  the  fundamental  frequency  of  vibration 
after  the  integration  is  completed. 
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2.4  Theoretical  Results 

The  natural  frequencies  of  vibration  are  expressed  in 
the  form  o>/p,  where  cd  is  the  natural  frequency  of  vibration 
and  p  is  the  fundamental  frequency  of  vibration  of  the  cor¬ 
responding  prismatic  beam,  given  by 


(2.2.5) 


P 


2.4-1  Convergence  of  the  Results 

To  illustrate  how  the  results  obtained  using  Galerkinrs 
method  converge,  the  fundamental  frequency  of  a  linearly 
tapered  beam  was  calculated  using  the  one,  two,  and  three 
term  trial  solutions,  and  the  frequency  of  the  second  mode 
was  calculated  using  the  two  and  three  term  solutions.  The 
values  obtained  for  each  mode  were  then  compared  to  each 
other,  and  to  the  values  obtained  from  the  closed  form  solu¬ 
tions.  Table  2-1  gives  the  values  obtained  for  the  funda¬ 
mental  frequency  of  vibration  of  linearly  tapered  beams  when 
the  one,  two  and  three  term  trial  solutions  were  used  in  the 
application  of  Galerkin ' s  method  to  the  equation  of  motion. 
These  results  are  presented  graphically  in  Figure  2-3. 

Table  2-2  and  Figure  2-4  give  the  values  obtained  for  the 
frequency  of  the  second  mode  of  vibration  when  the  two  and 
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TABLE  2-1  FUNDAMENTAL  MODE  FREQUENCY  RATIOS 
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TABLE  2-2  SECOND  MODE  FREQUENCY  RATIOS 
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Figure  2"3  Convergence  of  Galerkin's  Method  to  the  Fundamental  Frequency  of  Vibration 
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P 


Figure  2~4  Convergence  of  Gaierkin's  Method  to  the  Second  Natural  Frequency  of  Vibration 
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three  term  trial  solutions  were  used  in  Galerkin ' s  method. 

In  Table  2-3,  the  values  obtained  for  the  natural 
frequencies  of  vibration  when  the  one,  two,  and  three  term 
trial  solutions  were  used  in  the  application  of  Galerkin '  s 
method  to  the  equation  of  motion  for  a  linearly  tapered 
beam  are  compared  to  the  values  obtained  from  the  closed 


form  solutions  of  Timoshenko 


,  Ward 


,  and  Wrinch 


From  these  results,  the  following  conclusions  can  be 
made  about  the  accuracy  of  the  values  which  can  be  obtained 
by  using  this  method.  For  beams  of  constant  thickness 
0=0)  and  linearly  varying  width,  Galerkin' s  method  gives 
the  natural  frequency  of  the  2nd  mode  to  within  0.1  per  cent 
of  the  value  obtained  from  the  closed  form  solution  when 
three  terms  are  used  in  the  trial  solution.  For  beams 
with  tapers  in  both  directions,  acceptable  values  for  the 
natural  frequency  of  the  second  mode  can  be  calculated  for 
all  values  of  a  when  (3  is  less  than  0.8,  if  the  three  term 
trial  solution  is  used.  For  values  of  (3  greater  than  0.8 
the  slope  of  the  curves  of  co/p  vs  (3  (Figure  2-2)  changes 
from  negative  to  positive,  and  the  error  in  the  calculated 
value  of  the  frequency  increases. 


2.4-2  Natural  Frequencies  of  Vibration  of  Linearly  Tapered 
Beams 


The  values  of  the  frequencies  of  vibration  of  the 
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■  >  .ir  -  1  >  1 

991(W 


TABLE  2-3  NATURAL  FREQUENCIES  OF  VIBRATION  FOR  THE  LIMITING  CASES 
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first  three  modes  were  calculated  using  the  three  term 
trial  solution „  These  results  are  presented  in  tabular 
form  in  Table  2-4,  and  graphically  in  Figures  2-5,  2-6, 
and  2-7  *  The  broken  sections  of  the  curves  in  Figures 
2-6  and  2-7  are  for  the  values  of  (3  for  which  the  calcu¬ 
lated  frequencies  become  less  accurate. 

From  Figure  2-6  it  can  be  observed  that  for  any  fixed 
value  of  alpha  the  second  mode  frequency  varies  as  1/P,  for 
values  of  (3  less  than  0.5.  If  these  linear  sections  are 
extrapolated  to  the  values  of  (3  =  1.0,  for  the  limiting 
cases  of  a  =  0.0  and  a  =  1.0  the  values  obtained  for  cn/p 
are  within  five  per  cent  of  the  values  obtained  from  the 
closed  form  solutions.  When  a  similar  procedure  is  used 
to  extrapolate  the  values  of  the  third  mode  frequencies  from 
(3  =  0.2,  the  frequencies  obtained  for  (3  =  1.0  compare  much 
better  with  the  corresponding  values  obtained  from  the 
closed  form  solutions  than  the  values  obtained  by  using 
only  Galerkin 1 s  method. 

2.4-3  Natural  Frequencies  of  Vibration  When  Shearing  Force 
and  Rotatory  Inertia  are  Considered 

When  the  effects  of  shearing  force  and  rotatory  inertia 
are  considered,  the  common  factor  pA^t^/EIg  cannot  be  taken 
out  of  the  equation  of  motion,  and  hence  the  values  of  au/p 
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Figure  2“5  Fundamental  Frequency  of  Vibration  Versus  Thickness  Toper 
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Figure  2"6  Second  Natural  Frequency  of  Vibration  Versus  Thickness  Taper 
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Figure  2“  7  Third  Natural  Frequency  of  Vibration  VersusThickness  Taper 
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vary  with  the  mechanical  properties  of  the  material  used 
in  the  beam.  The  results  presented  in  this  section  are  for 

it 

beams  which  are  made  out  of  Alclad  7075-T6  aluminum  .  The 
mechanical  properties  of  this  material  are i 

Young's  Modulus  E  =  10.4  x  106  lb/in2. 

Shear  Modulus  G  =  3.9  x  10^  lb/in2. 

3 

Weight  Density  pg  =  0.101  lb/in  . 

The  dimensions  of  the  beams  are  taken  to  be  three 
inches  wide  and  one  quarter  inch  thick  at  the  clamped  end, 
and  the  lengths  considered  vary  from  two  inches  to  16  inches. 
When  the  effects  of  shearing  force  and  rotatory  inertia  are 
considered  in  a  prismatic  beam,  Equation  (2.2.9)  has  to  be 
solved  to  determine  the  natural  frequencies  of  vibration. 

The  natural  frequencies  of  the  first  three  modes  were  de¬ 
termined  from  this  equation,  and  are  compared  to  the  natural 
frequencies  which  were  obtained  when  shearing  force  and 
rotatory  inertia  were  neglected  in  Table  2-5. 

From  Table  2-5,  it  can  be  seen  that  shearing  force  and 
rotatory  inertia  have  little  effect  on  the  fundamental  fre¬ 
quency  of  vibration  of  beams  of  this  cross  section.  For 
the  second  mode,  the  difference  between  the  natural  frequencies 
of  vibration  which  were  calculated  when  shearing  force  and 


*  This  is  the  material  which  was  used  to  make  the  test 
specimens . 
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TABLE  2-5  FREQUENCIES  QF  VIBRATION  OF  UNIFORM  BEAMS 


LENGTI 

INCHES 

FREQUENCY  CPS . 

FIRST  MODE 

SECOND  MODE 

THIRD 

MODE 

SIMPLE 

SH.F . AND 

SIMPLE 

SH.F. AND 

SIMPLE 

SH.F. AND 

THEORY 

ROT . INR . 

THEORY 

ROT . INR . 

THEORY 

ROT . INR . 

16 

31.45 

31.45 

197.1 

197.0 

551.9 

550.8 

14 

41.08 

41.08 

257.4 

257.3 

720.9 

719.0 

12 

55.91 

55.91 

350.4 

350.0 

981.2 

977.8 

10 

80.51 

rH 

LO 

• 

o 

CD 

504.7 

503.8 

1412.9 

1405.8 

8 

125.8 

125.8 

788.4 

786.6 

2207.7 

2190.5 

6 

223.6 

223.6 

1402 

1398 

3925 

3871 

4 

503.2 

503.2 

3154 

3127 

8831 

8570 

2 

2012 

2012 

12608 

12189 

35326 

3.1786 
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rotatory  inertia  were  and  were  not  considered  varies  as 
approximately  1/1,  to  a  maximum  difference  of  -3.3  per 
cent  of  the  frequency  given  by  the  simpler  theory,  for  a 
two  inch  beaitio  Shearing  force  and  rotatory  inertia  have 
a  more  pronounced  effect  on  the  natural  frequency  of  the 
third  mode.  These  effects  reduce  the  frequency  from  0„2 
per  cent  for  a  16  inch  beam  to  10  per  cent  for  a  two  inch 
beam. 

The  fundamental  frequency  of  vibration  of  wedges  of 
lengths  between  two  and  16  inches  and  tapers  varying  from 
P  =  0o0  to  P  s  1,0  was  calculated  by  applying  the  one  term 
trial  solution  to  Equation  (2.3.13)  (See  Appendix  i-B) . 

The  results  which  were  obtained  are  presented  in  the  form 
of  frequency  ratios  (o>/p)  in  Figure  2“-8J  and  Table  2-6 . 

The  frequency  ratios  which  were  obtained  for  the  wedges 
between  eight  and  16  inches  long  were  all  equal  (for  equal 
tapers) ,  and  lower  than  the  corresponding  frequency  ratios 
obtained  when  the  three  term  trial  solution  was  applied  to 
the  simpler  equation  where  shearing  force  and  rotatory  in¬ 
ertia  are  neglected.  These  frequency  ratios  are  given  by  the 
solid  line  in  Figure  2-8 »  For  wedges  between  eight  and  two 
inches  long,  the  frequency  ratios  varied  slightly,  being 
higher  than  the  frequency  ratios  for  the  longer  wedges  in 
the  regions  near  {3  «  0.0  and  p  35  1.0,  and  lower  than  the 
frequency  ratios  for  the  longer  wedges  for  values  of  (3  between 
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0.06  and  0.94.  The  frequency  ratios  for  the  two  inch 
wedges  are  given  by  the  broken  line  in  Figure  2-8,  and 
the  frequency  ratios  for  the  wedges  between  two  and  eight 
inches  long  are  all  in  the  shaded  area  between  the  solid 
and  broken  lines. 

The  frequency  ratios  which  were  obtained,  for  the 
wedges  between  eight  and  16  inches  long,  when  the  effects 
of  shearing  force  and  rotatory  inertia  were  considered  are 
compared  to  the  values  obtained  from  the  simpler  case  when 
shearing  force  and  rotatory  inertia  were  neglected  in 
Table  2-6 „ 

TABLE  2-6  FUNDAMENTAL  FREQUENCY  RATIO  FOR  WEDGES 
BETWEEN  EIGHT  AND  16  INCHES  LONG 


BETA 

0.0 

0.2 

0.4 

0.6 

00 

0 

o 

1.0 

ONE 

TERM 

1  .  004 

1.043 

1.105 

1.206 

1.386 

1.785 

THREE 

TERMS 

1.000 

1.026 

1.063 

1.120 

1.226 

1.532 

SH «,  F. 

R  J, 

1.004 

1.026 

1.045 

1.064 

1 . 080 

1.096 

From  the  values  given  in  Table  2-6,  it  can  be  seen  that  the 
effects  of  shearing  force  and  rotatory  inertia  partially 
counteract  the  effects  of  taper  on  the  fundamental  frequency 


of  vibration. 
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Figure  2“8  Fundamental  Frequency  of  Vibration  of  Aluminum  Wedges  When  Shearing  Force 
and  Rotatory  Inertia  are  Conidered 
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2 .4-4  Effects  of  Nonlinear  Tapers 

The  one  term  trial  solution  was  used  to  determine 
what  effect  nonlinear  tapers  have  on  the  fundamental  fre¬ 
quency  of  vibration  of  cantilever  beams.  The  results  of 
this  problem,,  and  the  effects  of  nonlinear  tapers  on  the 
shapes  of  the  beams  are  presented  in  Figures  2-9  and  2-10 „ 

The  curves  in  Figure  2-9  are  an  indication  of  the  ef¬ 
fects  of  nonlinear  tapers  on  the  fundamental  frequency  of 
vibration.  The  values  obtained  for  the  case  of  n  -  2.G, 
m  —  2.0  have  approximately  the  same  accuracy  as  those  ob¬ 
tained  for  the  linearly  tapered  beams  (n  -  1.0,  m  ~  1.0), 
when  the  one  term  trial  solution  was  used.  For  the  case 
of  n  ”  0.5,  m  ™  0.5  the  abnormal  appearances  of  the  curves 
cast  some  doubt  on  the  reliability  of  the  results.  From 
Figure  2=10  it  can  be  seen  that  the  surfaces  of  this  beam 
are  concave,  and  have  a  very  large  curvature  near  the  root 
of  beam.  Because  of  this  it  becomes  difficult  to  predict 
the  behavior  of  the  beam. 

For  the  values  of  n  and  m  greater  than  one,  the  beam, 
approaches  the  shape  of  a  uniform  beam,  and  the  trial  solu¬ 
tion  is  a  more  accurate  representation  of  the  actual  deflected 
form.  This  should  give  more  accurate  results.  It  can  also 
be  expected  that  the  beam  will  behave  more  like  a  uniform 
beam  in  the  higher  modes  of  vibration. 
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Figure  2“9  Fundamental  Frequency  of  Vibration  Versus  Thickness  Taper  for  Nonlinearly 
Tapered  Beams 
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Figure  2_I0  Beam  Profiles 
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2  o  5  Summary 

In  this  Chapter,  Galerkin 0 s  method  was  used  to  calcu¬ 
late  the  first  three  natural  frequencies  of  lateral  vibration 
of  linearly  tapered  cantilever  beams.  The  frequencies  ob- 
tained  were  compared  to  the  frequencies  obtained  from  the 
closed  form  solutions  for  the  special  cases  of  a  uniform 
beam,  and  beams  which  have  an  area  of  zero  at  the  free  end. 

It  was  found  that  Galerkin ' s  method  converged  very  fast  to 
the  fundamental  frequencies  of  vibration,  and  when  the  three 
term  trial  solution  was  used,  the  value  calculated  for  the 
fundamental  frequency  of  a  pointed  beam  was  less  than  1.25 
per  cent  higher  than  the  frequency  given  by  the  closed  form 
solution  for  the  same  case. 

The  values  obtained  for  the  second  natural  frequencies, 
when  the  three  term  trial  solution  was  used,  were  much  better 
than  the  corresponding  values  obtained  when  the  two  term 
trial  solution  was  used,  but  for  values  of  (3  greater  than 
0o8s  the  errors  were  still  quite  large,  i.e.  30  -  40  percent. 

Since  at  least  n  terms  are  required  in  the  trial  solu- 

ttl 

tion  in  order  to  compute  a  value  for  the  n  natural  fre¬ 
quency,  only  a  first  approximation  was  obtained  for  the 
third  natural  frequencies,  and  the  frequencies  which  were 
calculated  were  much  too  high  to  be  of  any  practical  use. 
Better  values  could  be  obtained  for  the  third  natural 
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frequency  if  more  than  three  terms  were  used  in  the  trial 
solution » 

When  the  three  term  trial  solution  was  used,  the  de¬ 
terminant  of  a  three  by  three  symmetric  matrix  had  to  be 
set  equal  to  zero  to  determine  the  natural  frequencies  of 
vibration.  When  the  determinant  was  expanded,  it  became 
necessary  to  find  the  difference  between  large  numbers,  and 
the  round  off  errors  which  were  introduced  were  sufficient 
to  give  erroneous  results.  To  avoid  this,  the  matrix  was 
changed  by  means  of  a  series  of  matrix  multiplications 
and  the  frequencies  were  found  by  using  an  iterative  procedure. 

Also  considered  in  this  Chapter  were  the  problems  of 
finding  the  natural  frequencies  of  vibration  of  linearly 
tapered  cantilever  beams  when  shearing  force  and  rotatory 
inertia  were  considered,  and  the  natural  frequencies  of 
vibration  on  nonlinearly  tapered  beams. 

For  the  case  when  shearing  force  and  rotatory  inertia 
were  considered,  it  was  found  that  the  two  equations  which 
govern  the  motion  could  be  reduced  to  one  equation  for  beams 
with  constant  width  and  linearly  varying  thickness.  When 
Galerkin 1 s  method  was  applied  to  this  equation  using  the 
one  term  trial  solution,  the  values  obtained  for  the  funda- 
mental  frequency  were  much  lower  than  the  values  obtained 
when  the  effects  of  shearing  force  and  rotatory  inertia 
were  neglected.  In  view  of  this  fact  the  results  have  to 
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be  assumed  to  be  inaccurate,  however,  the  results  did  in¬ 
dicate  that  tapering  the  beams  amplify  the  effects  of  shear¬ 
ing  force  and  rotatory  inertia . 

Some  experiments  were  conducted  with  linearly  tapered 
cantilever  beams  to  determine  experimental  values  for  the 
first  three  natural  frequencies  of  vibration..  In  Chapter 
IV  these  experiments  are  discussed  and  the  experimental 
results  are  presented  and  compared  to  the  theoretical  results . 
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CHAPTER  III 


EXPERIMENTS  AND  EXPERIMENTAL  RESULTS 


Frequency  response  spectrums  were  obtained  for  eight 
and  16  inch  long  aluminum  beams  by  mounting  each  beam  on 
an  electrodynamic  shaker  table  and  recording  the  amplitudes 
of  the  deflections  of  the  tip  and  root  of  the  beam  at 
various  frequencies,, 

3 . 1  Apparatus 

The  vibrations  in  the  beams  were  excited  by  an  Unholtz- 
Dickie  electrodynamic  shaker ,  which  was  controlled  by  an 
automatic  vibration  exciter  control.  With  this  system, 
the  frequency  of  the  vibration  and  the  amplitude  of  the 
support  deflection  could  be  set  at  any  desired  level  which 
was  within  the  limits  of  the  testing  apparatus.  A  frequency 
meter  was  connected  to  automatic  vibration  exciter  control 
and  the  frequency  of  the  vibration  was  read  from  it.  The 
amplitude  of  the  support  deflection  was  read  directly  from 
the  meter  of  the  vibration  exciter. 

To  measure  the  tip  deflection,  a  stroboscope  and  a 
Gaertner  cathetometer  were  used.  The  stroboscope  had  to 
be  synchronized  with  the  frequency  of  the  shaker  and  the 
phase  angle  between  the  input  signal  and  the  output  of  the 
stroboscope  had  to  be  varied  through  360  degrees  in  order 
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to  determine  the  extreme  positions  of  the  tip  of  the  beam. 

The  frequency  synchronizing  and  phase  shifting  was  done 
automatically  with  a  Chadwick~Helmuth  Slip-Sync  system. 

A  signal  from  the  vibration  exciter  was  fed  to  the  Slip- 
Sync  system,  which  produced  a  pulse  output  at  a  frequency 
proportional  to  the  input  frequency  and  out  of  phase  with 
the  input  signal  by  an  angle  which  could  be  varied  from 
zero  to  360  degrees „  The  output  from  the  Slip-Sync  system 
was  used  to  control  the  stroboscope.  The  Slip-Sync  system 
was  set  to  vary  the  phase  angle  continuously  through  360 
degrees,  and  when  the  tip  of  the  beam  was  observed  in  the 
light  from  the  stroboscope,  it  appeared  to  be  oscillating 
slowly,  and  the  amplitudes  of  the  oscillations  were  measured 
with  the  cathetometer . 

3„2  Test  Specimens 

The  beams  were  machined  from  a  sheet  of  Alclad  7075-T6 
aluminum.  The  dimensions  of  the  root  were  three  inches 
wide  by  one  quarter  inch  thick,  and  the  lengths  were  eight 
inches  and  16  inches.  To  maintain  the  boundary  conditions 
at  the  clamped  end  of  the  beam,  (y(0)  =  0  and  *  0) 

the  beams  were  made  symmetric  about  the  clamped  area. 

This  produced  two  beams  of  identical  geometry  (See  Figure  3-1). 
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Figure  3 “I  Specimen  Configuration  and  Mounting 
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3„3  Experimental  Procedure 

The  frequency  response  spectrums  for  the  first  three 
modes  of  vibration  were  expected  to  look  like  the  frequency 
response  spectrum  in  Figure  3- 2, 


P 

FIG0  3“2  FREQUENCY  RESPONSE  SPECTRUM 

where  y/Yq  is  the  ratio  of  the  tip  deflection  to  the  sup¬ 
port  deflection called  the  magnification  factor ,  and  the 
values  of  co/p  which  correspond  to  the  peaks  A,  B,  and  C 
give  the  natural  frequencies  of  the  first  three  modes  of 
vibration o 

The  first  step  of  each  test  was  to  determine  the  lo¬ 
cation  of  the  first  three  peaks .  The  frequency  was  in¬ 
creased  slowly  from  zero  cycles  per  second  to  the  frequency 
at  which  the  third  peak  occurred,  and  the  frequencies  at 
which  the  peaks  occurred  were  recorded.  Following  this  step, 
the  tip  and  root  deflections  were  measured  at  various 
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frequencies  between  15  cps  and  1100  cps  for  the  eight  inch 
beams,  and  15  cps  and  500  cps  for  the  16  inch  beams.  It 
was  initially  planned  to  measure  the  tip  deflections  at 
all  frequencies  for  a  constant  support,  deflection  of  0.02 
inches,  but  some  difficulties  were  encountered,  and  the 
support  deflections  had  to  be  changed  during  each  test,. 

At  frequencies  near  the  natural  frequencies  the  tip 
deflections  were  very  large,  and  the  amount  of  time  that 
the  tip  of  the  beam  spent  in  the  vicinity  of  the  extreme 
points  of  its  movement  was  very  small,  and  it  was  difficult 
to  locate  maximum  and  minimum  positions  for  the  tip  of  the 
beam a  To  reduce  the  tip  movement  the  support  deflection 
was  reduced  to  0.QQ1  inches,  and  more  accurate  measurements 
of  the  tip  movement  could  be  made*  Wien  using  this  very 
small  support  deflection  it  became  difficult  to  measure 
the  tip  deflection  at  frequencies  which  were  not  near  the 
natural  frequencies 0  At  these  frequencies,  and  with  a  sup- 
port  deflection  of  0.001  inches,  the  tip  of  the  beam  was 
moving  less  than  0.QQ4  inches.  This  small  movement  was 
hard  to  detect  and  could  not  be  measured  accurately,  there¬ 
fore  larger  support  deflections  had  to  be  used  for  frequencies 
which  were  in  this  range. 

Another  problem  arose  from  the  fact  that  the  force 
required  to  vibrate  the  beam  varies  as  the  product  of  the 
support  deflection  and  the  square  of  the  frequency.  Because 
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of  this  fact,  the  force  required  to  vibrate  the  beam  at 
high  frequencies  was  greater  than  the  shaker  could  produce „ 

To  expand  the  range  of  frequencies  over  which  readings  could 
be  taken,  the  support  deflections  were  decreased  as  the 
frequency  was  increased,  until  the  tip  movements  became  too 
small  to  be  measured  with  the  existing  apparatus „  As  a 
result.,  the  frequency  spectrums  could  not  be  completed  for 
the  third  modes  of  the  eight  inch  beams „  However,  the  nat¬ 
ural  frequencies  of  the  third  modes  could  still  be  determined 
because  the  tip  movements  became  large  enough  to  be  seen 
when  the  third  natural  frequency  was  reached,, 

When  it  became  necessary  to  vary  the  support  deflections, 
a  series  of  tests  were  conducted  to  determine  whether  or  not 
the  magnification  factors  varied  linearly  with  the  support 
deflections „  These  tests  and  the  results  obtained  from  them, 
are  discussed  in  Subsection  303-l0 

Since  the  support  deflections  had  to  be  varied,  the 
tests  were  conducted  in  the  following  manner 0  The  initial 
frequency  for  each  test  was  15  cps,  and  the  frequency  was 
increased  by  amounts  varying  from  two  cps  for  frequencies 
which  were  near  the  natural  frequencies  to  50  cps  for  fre¬ 
quencies  which  were  far  from  the  natural  frequencies  (i0e0  in 
the  regions  where  the  slope  of  the  frequency  spectrum  was 
near  zero) „  The  support  and  tip  deflections  were  measured 
at  each  frequency  when  the  support  deflection  was  such  that 
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the  tip  deflection  was  in  the  range  where  it  could  he  mea¬ 
sured  most  accurately  (0.01  to  0.4  inches) .  The  frequency 
spectrums  obtained  in  this  manner  are  presented  in  Section 
3  o  4. 

The  values  of  the  natural  frequencies  obtained  from 
these  tests  were  always  lower  than  the  theoretical  values 
(See  Table  3-2) ,  and  when  the  beams  were  examined  after  each 
test s  abraded  areas  were  found  in  the  clamped  area  of  the 
beams.  This  indicated  that  there  was  relative  movement  be¬ 
tween  the  clamp  and  the  beam,  which  could  occur  if  the  clamp 
was  loose  or  deforming  with  the  beams.  Tests  were  then 
conducted  to  determine  the  effect  of  clamping  pressure  on 
the  fundamental  frequency  of  vibration  of  the  beam  (See 
Subsection  3.3-2).  It  was  found  that  by  changing  the  clamp¬ 
ing  pressure,  no  appreciable  change  in  the  fundamental  fre- 
quency  could  be  induced.  The  possibility  of  the  clamp  de¬ 
forming  with  the  beam  is  discussed  in  Section  3-4. 

3.3-1  Linearity  of  the  System 

In  order  to  maintain  tip  deflections  which  could  be 

measured  accurately,  it  was  necessary  to  vary  the  support 

deflections  as  the  frequency  was  changed.  A  series  of  tests 

was  conducted  to  determine  whether  or  not  the  magnification 

factor  )  remained  constant  as  the  support 

'root  deflection 
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deflection  was  varied  through  the  range  of  values  for  which 
the  tip  deflections  could  be  measured  at  a  constant  fre¬ 
quency.  A  typical  set  of  results  is  presented  in  Figure 
3“3=.  From  Figure  3-3,  it  can  be  seen  that  at  frequencies 
which  are  not  near  a  natural  frequency,  the  magnification 
factor  (y/Yq)  decreases  slightly  as  the  support  deflections 
are  increased.  For  the  cases  when  the  frequencies  were 
near  the  natural  frequencies  of  vibration  (in  this  case  30 
cps  for  the  first  mode  and  180  cps  for  the  second  mode) 
the  magnification  factors  were  large  and  unstable „  In  these 
regions  the  tip  deflections  would  vary  considerably  when 
the  support  deflection  and  frequency  were  left  constant, 
and  the  tip  deflections  could  not  be  measured  accurately » 

303“2  Effects  of  Clamping  Pressure  on  the  Fundamental 
Frequency 

The  abrasions  found  on  the  clamped  area  of  the  beams, 
indicated  that  the  beams  were  moving  in  the  clamp.  If  the 
movement  between  the  clamp  and  the  beam  could  be  reduced 
by  increasing  the  clamp  pressure,  the  natural  frequency  of 
vibration  of  the  beam  would  increase  because  the  effective 
length  of  the  beam  would  be  reduced.  A  series  of  experiments 
was  conducted  to  determine  what  effect  the  clamp  pressure 
had  on  the  natural  frequency  of  vibration. 
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Figure  3~3  Magnification  Factor  Versus  Support  Deflection  fora  16  Inch  Uniform  Beam 
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To  record  the  vibrations  of  the  beam,  an  electric 
resistance  strain  gauge  was  mounted  near  the  root  of  the 
beam,  and  connected  to  an  oscillograph „  The  beam  was  put 
in  the  clamp,  and  the  bolts  on  the  clamp  were  tightened 
with  a  torque  wrench.  The  free  end  of  the  beam  was  given 
a  small  deflection  and  released,  and  the  vibration  was  re¬ 
corded.  This  was  repeated  for  four  different  values  of 
bolt  torque  with  the  eight  and  16  inch  uniform  beams 0  The 
results  of  these  tests  are  given  in  Table  3-1 . 


TABLE  3-1  EFFECT  OF  CLAMP  PRESSURE  ON  THE 

FUNDAMENTAL  FREQUENCY  OF  UNIFORM  BEAMS 


BOLT  TORQUE  NATURAL  FREQUENCY  (CPS) 

(IN.  LB.)  8  INCH  16  INCH 


80 

100 

125 

150 


116.16 

30 . 14 

116.30 

30.16 

116.54 

30.17 

116.67 

30.17 

3.3-3  Discussion  of  the  Experiments 


The  natural  frequencies  of  vibration  were  assumed  to 
be  those  frequencies  at  which  the  magnification  factors 
reached  a  maximum.  This  assumption  is  correct  when  there 
is  no  damping,  but  when  damping  is  present,  the  frequency 
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at  which  the  maximum  magnification  occurs  is  o>/p  -  h  -  252  9 
and  the  natural  frequency  is  0)/p  ~  ’ll  -  6 2  ,  where  6  is 
the  damping  factor.  From  the  expressions  for  0)/p,  it  can 

be  seen  that  for  small  values  of  6  the  natural  frequency  is 
very  close  to  the  frequency  at  which  the  maximum  magnification 
occurs.  Since  the  damping  was  assumed  to  be  small  for  the 
beams  being  tested^ the  assumption  that  the  natural  frequency 
is  the  frequency  at  which  the  maximum  magnification  occurs 
is  valid. 

In  view  of  the  results  of  the  tests  conducted  for 
Subsection  3„3-ls  it  must  be  kept  in  mind  that  when  viewing 
the  frequency  spec t rums  the  values  of  y/yQ  greater  than 
seven  are  possibly  in  error ,  because  of  the  instability  of 
the  system  when  the  magnification  factors  are  large.  How- 
ever,  this  does  not  have  a  serious  effect  on  the  results 
which  are  desired  from  the  frequency  spectrums,  since  the 
main  object  of  these  experiments  was  to  determine  the  natural 
frequency  of  vibration ,  and  not  the  maximum  magnification 
factors  which  could  be  obtained. 

3.4  Experimental  Results 

The  frequency  speetrums  presented  in  this  section  were 
obtained  by  measuring  the  amplitudes  of  the  tip  deflections 
(y)  when  the  vibrations  were  excited  by  a  known  support 
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deflection  (yQ) .  The  magnification  factors  are  defined 
as  Y/Yq  an<^  the  frequencies  are  given  in  the  form  cu/p, 
where  cd  is  the  exciting  frequency  and  p  is  the  fundamental 
frequency  of  the  corresponding  uniform  beam.  The  frequency 
spectrums  obtained  in  this  manner  are  presented  in  Figures 
3- 4  to  3-1 03  and  the  experimental  values  of  the  natural 
frequencies  are  given  in  Table  3-2, 

3,4-1  Discussion  of  the  Experimental  Results 

On  the  basis  of  the  experimental  results  which  were 
obtained  and  observations  which  were  made  while  the  experi¬ 
ments  were  being  conducted,  the  following  conclusions  were 
made  about  the  experimental  results. 

The  natural  frequencies  obtained  experimentally  were 
lower  than  the  values  given  by  Galerkin’s  method,  and  also 
generally  lower  than  the  values  which  were  computed  using 
Martin’s  method.  It  is  also  apparent  that  the  longer  beams 
gave  better  experimental  values  than  the  shorter  beams. 

On  examining  each  beam  after  the  tests  were  completed,  an 
abraded  area  was  found  near  the  edge  of  the  clamp.  This 
indicated  that  there  was  relative  movement  between  the  sur¬ 
faces  of  the  clamp  and  the  surfaces  of  the  beam.  Since  the 
clamp  pressure  did  not  affect  the  natural  frequency  of  vibra¬ 
tion  (See  Subsection  3,3-2),  the  clamp  had  to  be  deforming 
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Figure  3“4  Frequency^ Spectrum  for  an  8" Aluminum  Beam  c*=0.0,  p=00,  p=  125.9  cps 
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Figure  3~5  Frequency  Spectrum  fora  !6"Aluminum Beam  a=0j0,p=0.0>p=  31.5 cps 
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Figure  3-6  Frequency  Spectrum  for  an  8"  Aluminum  Beam  a=0.0, 8=0.8,  p=  125.9  cps 
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Figure  3“  7  Frequency  Spectrum  for  a  16" Aluminum  Beam  oi=0.0,  p=0.8,  p  =  3l.5  cps 
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Figure  3~8  Frequency  Spectrum  for  an  8“Aluminum  Beam  a=0.8,  p=00,p=l25.9cps 
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TABLE  3-2  EXPERIMENTAL  VALUES  FOR 
FREQUENCY  RATIOS  -  a>/p 


METHOD 

LENGTH 

TAPER 

MODE 

GALERKIN ' S 

MARTIN’S* 

EXPERIMENTAL 

CLOSED 

INCHES 

a 

3  TERMS 

FORM 

16 

o 

0 

° 

o 

c 

o 

1 

1.000 

1.000 

0.95 

1.000 

p=31.5  cps 

2 

6.268 

6.267 

5.90 

6.267 

3 

18.846 

17.551 

16.4 

17.551 

8 

o 

0 

o 

o 

o 

o 

1 

1.000 

1.000 

0.93 

1.000 

pss125.9  cps 

2 

6.268 

6.267 

5.5 

6.267 

3 

18.846 

17.551 

15.5 

17.551 

16 

o 

e 

00 

0.0 

1 

1.535 

1.368 

1.5 

p=31.5  cps 

2 

7.307 

6.994 

6.8 

3 

21.155 

18.222 

17.5 

8 

00 

0 

o 

o 

0 

o 

1 

1.535 

1.368 

1.4 

p-125 . 9  cps 

2 

7.307 

6.944 

6.7 

3 

21.155 

18.222 

16.6 

16 

o 

0 

© 

00 

c 

o 

1 

1.226 

1.140 

1.2 

p-31.5  cps 

2 

4.897 

4.538 

4.8 

3 

17.615 

10.892 

11.5 

8 

0.0 

0.8 

1 

1.226 

1.140 

1.1 

p-125.9  cps 

2 

4.897 

4.538 

3.8 

3 

17.615 

10.892 

8.7 

16 

0.5 

0.5 

1 

1.316 

1.277 

1.2 

p-31.5  cps 

2 

5.693 

5.551 

4.7 

3 

18.788 

13.865 

11.9 

* 


The  values  calculated  by  Martin 


give  a  lower  bound 


to  the  theoretical  frequencies  for  tapered  beams. 
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with  the  vibrating  beam.  With  the  clamp  deforming,  the 
point  where  the  boundary  conditions,  y  *=  0  and  -  0, 
were  satisfied  had  to  move  to  a  distance  greater  than  t 
from  the  free  tip  of  the  beam.  For  the  case  of  an  eight 
inch  uniform  beam,  it  was  found  that  the  experimental  value 
for  the  fundamental  frequency  was  (ju/p  ~  0.93,  where  p  is 
given  by 


EI0  /.1,8/sV2 
PA0  \  1  J  ° 


(2.2.5) 


If  all  the  parameters,  except  h,  in  the  formula  for  p  are 
assumed  to  be  exact,  the  error  in  p  due  to  an  error  in  the 
length  is  given  by 

dp  _  2d-l 
P  h 

In  order  for  this  to  give  a  fundamental  frequency  ratio  of 
0.93  instead  of  1.00,  dt  must  be  equal  to  +0.28  inches. 

On  the  basis  of  the  marks  observed  on  the  beam,  a  dh  of 
0.28  inches  was  indeed  possible,  and  it  was  concluded  that 
the  deformation  of  the  clamp  was  an  important  factor  which 
could  cause  an  error  in  the  experimental  frequency. 

Another  observation  which  was  made  is  that  the  error 
between  the  experimental  values  and  the  theoretical  values 
increased  as  the  mode  of  vibration  increased,  and  also  that 
the  range  of  frequencies  around  each  mode  for  which  the 
magnification  factor  was  large  increased  for  each  higher 
mode  (See  Table  3-3) .  Since  the  amount  of  energy  in  the 
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system  increased  as  the  mode  of  vibration  was  increased, 
the  amount  of  energy  transferred  to  the  clamp  also  had  to 
increase „  This  would  give  larger  clamp  deflections  and • 
thereby  increase  the  effective  length  of  the  beam.  A  second 
factor  which  contributed  to  the  larger  error  at  the  higher 
modes  was  the  increased  effect  of  shearing  force  and  rota¬ 
tory  inertia,  which  also  tend  to  lower  the  natural  frequency 
of  vibration. 

From  the  fact  that  the  range  of  frequencies  for  which 
the  magnification  factor  was  large  increased  as  the  mode  of 
vibration  increased,  it  is  suspected  that  the  damping  present 
in  the  beams  also  increased  with  frequency.  Since  the  effect 
of  damping  is  to  decrease  the  natural  frequencies  of  vibra¬ 
tion,  the  increases  in  the  damping  may  also  have  contributed 
to  the  error  that  increased  with  the  mode  of  vibration. 

The  next  observation  made  was  that  for  two  tapered 
beams  made  to  the  same  specifications,  (end  A  and  end  B  of 
the  existing  beams)  the  natural  frequencies  of  vibration 
of  the  second  and  third  modes  varied  as  much  as  eight  per 
cent  between  the  two  beams,  with  the  largest  differences 
being  recorded  for  the  beams  with  a  thickness  taper  only. 

A  careful  examination  of  the  beams  revealed  that  the  top  and 
bottom  surfaces  were  not  plane.  Since  the  thickness  dimen¬ 
sion  of  the  beams  was  small  (0.25  inches  at  the  root  and  de¬ 
creasing  toward  the  free  end) ,  small  variations  in  the 
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thickness  dimension  could  have  had  a  significant  effect  on 
the  stiffness  of  the  beams  and  ultimately  effect  the  natural 
frequencies  of  vibration.  The  irregularity  in  the  surfaces 
of  the  beams  could  possibly  have  contributed  to  the  dif¬ 
ferences  in  the  natural  frequencies  of  two  beams  which  were 
made  to  the  same  specifications. 

The  next  observation  made  was  that  for  magnification 
factors  of  greater  than  y/y0  “  7.0,  the  tip  deflections  were 
unstable.  This  was  indicated  by  the  scatter  of  the  points 
plotted  on  the  frequency  spectrums  where  the  magnification 
factors  are  greater  than  7,0,  When  the  beams  were  vibrating 
at  a  frequency  near  a  natural  frequency,  the  variations  in 
the  tip  deflections  were  plainly  visible. 

In  Table  3-3  the  range  of  frequencies  around  each  mode 
for  which  the  magnification  factors  were  greater  than  6.0 
are  given.  This  table  illustrates  the  fact  that  the  range 
of  frequencies  around  each  mode  for  which  the  magnification 
factors  were  large  increased  as  the  mode  of  vibration  in¬ 
creased.  It  can  also  be  seen  from  this  table  that  the  fre¬ 
quencies  which  were  calculated  for  the  first  and  second  modes 
using  the  three  term  trial  solution  in  Galerkin“s  method  are 
usually  in  the  range  where  the  experimentally  obtained  mag¬ 
nification  factors  are  greater  than  6,0,  The  values  which 
were  computed  for  the  third  mode  are  a  first  approximation 
to  the  third  natural  frequency,  and  could  be  improved  if 
more  terms  were  taken  in  the  trial  solution. 
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TABLE  3-3  RANGE  OF  FREQUENCIES  FOR  WHICH  THE 

MAGNIFICATION  FACTORS  WERE  GREATER  THAN  6.0 


TAPER 

MODE 

RANGE  OF  co/p  FOR  y/y0 

>6.0 

CALCULATED 

a 

P 

1 

8  18  BEAM 

16 88 

BEAM 

a)/p 

1 

0 . 80  -  1.00 

0  o  75 

-  1.10 

loQQO 

cT 

0 

o 

o 

0 

o 

2 

5  o 09  -  6,01 

5  o  5 

-  6.3 

6.288 

3 

— 

15 .  7 

-  17.2 

18.846 

1 

1.20  -  1 .67 

1  „  2  - 

1.7 

1.535 

o 

0 

CD 

o 

0 

o 

2 

5  a  8  -  8 . 0 

6.1  - 

7.7 

7.307 

3 

— 

16.3  - 

18.9 

21.155 

1 

0.8  -  1.3 

0 . 95 

-  1.4 

1.226 

O 

0 

o 

CX) 

<> 

O 

2 

3.1  -  4  0  7 

4.1  - 

5.7 

4.897 

3 

— 

9  o  8  - 

13.4 

17.615 

1 

_ 

1.1  - 

1.4 

1.316 

0  .  5 

in 

o 

o 

2 

— 

4.2  ~ 

5.5 

5 . 693 

3 

— 

10.7  - 

13.3 

18.788 
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CHAPTER  IV 


CONCLUSIONS  AND  DISCUSSION 

4.1  Conclusions 

The  values  which  were  obtained  from  Galerkin 8 s  method 
for  the  natural  frequencies  of  vibration  are  an  upper  bound 
to  the  frequencies  which  were  obtained  from  the  closed  form 
solutions.  This  can  be  seen  in  the  results  which  are  pre- 
sented  in  Chapter  III. 

For  linearly  tapered  cantilever  beams ,  this  method  con¬ 
verged  very  well  to  the  fundamental  frequency  of  vibration, 
for  all  values  of  width  and  thickness  taper,  when  the  three 
term  trial  solution  was  used.  When  the  thickness  taper  was 
less  than  (3  -  0.5  the  two  term  trial  solution  gave  accept¬ 
able  values  for  the  fundamental  frequency.  The  three  term 
trial  solution  gave  acceptable  values  for  the  second  fre¬ 
quency  of  vibration  for  all  values  of  a  when  b  was  less 
than  0„8,  For  values  of  f3  greater  than  0.8  the  results  be¬ 
gan  to  diverge  and  the  error  increased  as  (3  approached  1.0. 

By  comparing  the  frequencies  obtained  from  Galerkin 8 s 

* 

method  to  the  frequencies  obtained  from  the  closed  form 
solutions,  it  was  found  that  the  values  obtained  for  the 
third  natural  frequency  were  7.7  per  cent  high  for  a  -  0.0 
and  |3  -  0.0,  and  11.6  per  cent  high  for  the  case  a  ~  1.0, 
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P  ®  0.0.  For  the  cases  where  (3  ~  1„0  the  frequencies 
obtained  were  over  IDO  per  cent  high0  These  errors  could 
be  reduced  if  more  terms  were  used  in  the  trial  solution's «, 
However  it  should  be  emphasized  that  a  degree  of  difficulty 
was  encountered  in  solving  the  determinant  which  was  ob- 
tained  when  the  three  term  trial  solution  was  used.  The 
difficulty  encountered  in  this  case  was  finding  the  dif¬ 
ferences  between  large  numbers  when  the  differences  were 
very  small .  This  condition  would  probably  get  worse  if 
more  terms  were  used  in  the  trial  solution.  With  this  fact 
in  mind,  it  was  concluded  that  Galerkin“s  method  is  not 
suitable  for  finding  the  frequencies  of  vibration  of  modes 
higher  than  the  third  mode  when  the  trial  solution  for  the 
deflected  shape  is  of  the  form 

y  *  a1(6X2-4X3+X4)  +  a ? ( 2QX2-X0X3+X5)  +  a3 (45X2-20X3+X6)  +  . . . 

From  the  values  obtained  for  the  second  natural  fre- 
quency  of  vibration  of  linearly  tapered  beams ,  it  was  noted 
that  for  beams  of  the  same  width  taper ,  the  decrease  in  the 
frequency  of  the  second  mode  was  very  nearly  linear  as  (3 
was  increased  from  0.0  to  0.5.  When  these  linear  sections 
were  extrapolated  to  j3  —  1.0  for  a  =  0.0  and  a  -  1.0,  the 
frequencies  obtained  for  the.  second  mode  agreed  very  well 
with  the  frequencies  which  were  obtained  from  the  closed  form 
solutions.  On  this  basis  it  was  concluded  that  the  frequency 
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of  the  second  mode  of  vibration  is  nearly  linearly  depend¬ 
ent  on  the  thickness  taper  on  the  beam.  This  is  discussed 
further  in  Section  4.3 „ 

It  was  found  that  the  effects  of  shearing  force  and 
rotatory  inertia  have  a  negligible  effect  on  the  fundamental 
frequency  of  vibrations  of  uniform  beams.  When  the  frequencies 
of  the  second  and  third  modes  were  calculated,  the  influence 
of  these  effects  was  found  to  increase  as  the  mode  of  vibra¬ 
tion  was  increased  and  as  the  length  of  the  beams  decreased. 
When  Galerkin ' s  method  was  applied  to  the  equation  of  motion 
for  a  wedge  when  the  effects  of  shearing  force  and  rotatory 
inertia  were  considered,  it  was  found  that  tapering  the  beam 
magnified  the  effects  of  shearing  force  and  rotatory  inertia 
on  the  fundamental  frequency  of  vibration. 

When  the  values  obtained  for  the  fundamental  frequency 
of  vibration  of  a  wedge  that  were  obtained  when  shearing 
force  and  rotatory  inertia  were  considered  were  compared  to 
the  frequencies  obtained  from  the  closed  form  solutions, 
and  the  experimental  frequencies,  it  was  found  that  for  a 
beam  with  a  -  0.0  and  £  ~  1.0,  the  frequency  obtained  was 
27.5  per  cent  lower  than  the  frequency  obtained  from  the 
closed  form  solution,  and  that  the  experimental  frequencies 
were  higher  than  the  theoretical  frequencies.  On  the  basis 
of  these  facts  it  was  concluded  that  the  values  obtained 
from  Galerkin1 s  method  for  this  case  were  in  error.  This  is 
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possibly  due  to  the  round  off  errors  that  were  present 
when  the  frequency  equation  was  being  solved. 

The  experimental  results  obtained  in  the  form  of  fre¬ 
quency  spectrums  show  that  the  range  of  frequencies  around 
each  mode  for  which  the  magnification  factors  were  large 
increased  as  the  mode  of  vibration  increased.  This  indicates 
that  the  damping  present  is  a  function  of  the  frequency. 

The  increasing  damping  partially  explains  the  larger  errors 
between  the  theoretical  results  and  the  experimental  results 
at  the  higher  modes  of  vibration. 

4.2  Summary  of  Conclusions 

1.  The  natural  frequencies  of  vibration  of  linearly  tapered 
beams  obtained  using  Galerkin's  method  are  higher  than  the 
experimental  frequencies  and  the  frequencies  obtained  from 
the  closed  form  solutions. 

2.  For  values  of  (3  less  than  0.5,  Galerkin's  method  gives 
acceptable  values  for  the  frequency  of  the  second  mode  when 
three  terms  are  used  in  the  trial  solution. 

3.  With  the  trial  solution  which  was  used,  Galerkin's 
method  is  not  suitable  for  calculating  the  natural  frequencies 
of  modes  higher  than  the  third  mode. 

4.  The  natural  frequencies  of  vibration  of  the  second  mode 
decrease  almost  linearly  as  the  thickness  taper  is  increased. 
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5.  When  the  effects  of  shearing  force  and  rotatory  in¬ 
ertia  are  considered,  Galerkin's  method  does  not  give 
satisfactory  values  for  the  fundamental  frequency  of  vibration, 

4,3  Second  Mode  Frequencies 

When  the  values  obtained  using  the  three  term  trial 
solution  were  compared  to  the  values  obtained  from  the  closed 
form  solutions ,  it  was  noted  that  the  frequencies  of  vibration 
of  the  second  mode  were  within  0ol  per  cent  for  beams  of 
constant  thickness  and  linearly  varying  width.  It  was  also 
found  that  for  beams  of  the  same  width  taper  the  second  nat¬ 
ural  frequency  decreased  almost  linearly  as  the  thickness 

•If 

taper  was  increased  from  p  =  CKO  to  p  =  0,5,  and  when  these 
linear  sections  were  extrapolated  to  p  =  1,0  for  a  -  Q.O  and 
ct  -  1,0  the  values  given  for  the  second  natural  frequencies 
were  within  five  per  cent  on  the  frequencies  given  by  the 
closed  form  solutions.  This  led  to  the  conclusion  that  the 
frequency  of  the  second  mode  of  vibration  decreases  linearly 

v 

as  the  thickness  taper  is  increased. 

When  the  second  natural  frequencies  which  were  computed 
from  the  closed  form  solutions  (See  Table  4-1)  were  compared, 
it  was  found  that  for  beams  with  a  ~  0,0  the  ratio  of  the 
frequency  when  P  -  1,0  to  the  frequency  when  p  -  0,0  was 
0.691,  and  for  beams  with  a  =  1,0  the  ratio  of  the  frequency 
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TABLE  4-1  SECOND  MODE  FREQUENCY  RATIOS  OBTAINED  FROM 


THE  CLOSED  FORM  SOLUTIONS 


0.0 


1.0 


CKO 


6.267 


4  0  3  26 


1.0 


8.830 


6.016 


when  (3  =  1.0  to  the  frequency  when  (3  =  0.0  was  0.681.  By 
assuming  that  the  values  of  the  corresponding  ratios  remain 
between  the  limits  0.681  and  0.691  for  all  values  of  a,  and 
using  the  fact  that  the  second  natural  frequency  decreases 
linearly  as  (3  increases,  it  has  been  found  that  the  frequency 
of  the  second  mode  of  vibration  can  be  estimated  accurately 
for  any  combination  of  linear  width  and  thickness  taper  once 
the  frequency  of  vibration  of  the  corresponding  beam  of  con¬ 
stant  thickness  has  been  determined. 

To  find  the  frequencies  it  was  assumed  that  for  any  fixed 
value  of  ctj  the  ratio  of  the  frequency  when  £  ~  1.0  to  the 
frequency  when  (3  -  0.0  was  0.686.  The  frequency  ratio  for 
any  value  of  a  and  (3  was  then  found  from  the  relationship 
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Figure  4”l  Second  Natural  Frequency  of  Vibration  Versus  ThicknessTaper 
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values  of  a  and  (3,  and  is  the  frequency  ratio  for  a 

beam  with  the  same  value  of  a  when  f3  =  0.0.  The  results 
obtained  from  this  method  are  superimposed  on  the  results 
obtained  when  Galerkin ' s  method  was  used  in  Figure  4-1, 
and  are  compared  to  the  experimentally  obtained  frequencies 
in  Table  4-2* 


TABLE  4“ 2  SECOND  MODE  FREQUENCY  RATIOS 


a 

p 

EXPERIMENTAL 

THEORETICAL 

16"  BEAM 

8"  BEAM 

LINEAR 

GALERKIN ' S 

o 

e 

o 

0.8 

4.8 

3.8 

4.68 

4.90 

0.8 

OoO 

6.8 

6.7 

7.31 

7.31 

o 

1  0 

U1 

0.5 

4.7 

- 

5.65 

5.69 

From  the  results  given  in  Table  4=2 ,  it  can  be  seen 
that  the  frequency  ratios  obtained  when  this  linear  extra¬ 
polation  was  used  are  closer  to  the  experimental  frequency 
ratios  than  the  frequency  ratios  obtained  by  using  Galerkin "s 
method  when  thickness  taper  was  considered 0  It  should  also 
be  noted  that  this  method  reduces  the  amount  of  work  required 
to  calculate  the  frequencies  of  the  second,  because  the 
equation  of  motion  reduces  from 
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i  t  i 

1  Y 


+ 


21 


9  19  1 

Y 


+  IY 


IV 


9  O 

pAY 


for  beams  with  tapers  in  both  directions  to 


E 


8  B 

21  Y 


IV 

+  IY 


oo 

35  -PAY 


for  beams  with  a  width  taper  only,  and  the  expressions  for 

/ 

the  area  and  moment  on  inertia  become  linear . 
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APPENDIX  I 

DERIVATION  OF  THE  EQUATION  OF  MOTION  FOR  A  WEDGE  WITH 
SHEARING  FORCE  AND  ROTATORY  INERTIA  CONSIDERED 0 


For  the  wedge  shown  in  Figure  1-1, 


t 

where  £3 

A 

and  I 


A0(l.  -  PX)  , 

I 0 ( 1  -  f3X)  3  ,  where  X  ~ 


When  shearing  force  and  rotatory  inertia  are  con¬ 
sidered  the  equations  which  govern  the  motion  ares 


d 

dx 


+  k  °  AG 


I  ip 


(2.2.9) 


k  AG 


A 

dx 


2  A 

-pco  Ay 


(2.2.10) 
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Expressing  A  as 


A 


and  defining 


D 


the  equations  become 


dx 


E(l-p.X) 


3  M 

dx 


0 


IQ  ( 1“>(3X) 


k  GA 


0 


b  ’ 


+  D(1=PX)(||  -  *)  -  -pcu2(l-pX)3  f  (1.1) 


d 

dx 


(1“PX) 


,  dx 


Let  (1-PX) 
then 


u 


__d_ 

dx 


(1~PX)  y 

k  G 


d  du 
du  dx 


(1*2) 


.  £  jL 

t  du 


u  du  e 

Substitute  u  into  equations  1.1  and  1.2  and  differentiate 
1.1  with  respect  to  u  to  get 


2n 
c  E 


s  2  11  11  ^  1  J  • 

6u ip  +  6u  ip  +  u  ip 


-  D 


du 


u  ( cy  +  ip) 


•pen 


2  3  11 

3  u  ip  +  u  ip 


d 

du 


u(cy  +  ip) 


pg> 


ck  G 


uy  , 


(1.3) 


(1*4) 


where  the  primes  represent  differentiation  with  respect  to  u3 


u 
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Substituting  the  right  hand  side  of  1.4  into  1.3  gives 

2 


2 

C  E 


1  2  0  0  O  0  0  0 

6uf  +  6u  ip  +  u 'ip 


Dpg) 

ck°G 


uy 


■pOi 


3  u2f  +  ip 


(1.5) 


Elimination  of  ip  from 


0  2  1  “  I!  I  0 

6uip  +  6u  ip  +  u  ip 


(A) 


i  e 


0  0  0  0  9 


du 


(1.4)  ”  2(cy  +  ip  )  +  u(ey  +  ip  ) 


^-7  (uy  +  y) 

ck  G 


2ip 


2cy  +  u(cy  +  ip  )  +  *^r  (uy  +  y) 

ck  G 


(I  c.6) 


du 


2  (1.4)  gives 


0  0  0 


Uf 


ooo  so  jy  nn\  8  8  8 

3  (cy  +  ip  )  +  ucy  +  (uy  +  2y  ) 

ck  G 


(1.7) 


Substitution  of  X„6  and  1.7  into (A) gives  a  term  which  is 
free  from  ipa  Equation  X.5  then. becomes 


c3E 


2  IV  ,  ,  '  '  '  ,  ^  '  '  ,  po) 

u  y  +  6uy  +  6y  +  — •%— 

-  “  “  c  k  G 


2  »  0  0 

u  y  +  5uy  +  3y 


^  2 
ck  G 


pen 


3uip  +  1 1  ip 


(1.8) 


Elimination  of  ip  from 


2  ' 

3uf  +  u  ip 


(B) 


. 
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Substituting  u  into  1.1  gives 


c2E 


i  2  1  1 

3u  ip  +  U  Ip 


•  2  2 

D  (cy  +  ip)  -  -po)  u  ip 


from  which  ip  can  be  found  in  terms  of  y  and  derivatives 
of  ip. 


2 

c  E 


ip  - 


8  2  1  1 
3u  ip  +  u  ip 


-  Dcy 


2  2x 

(D  -  pa)  u  ) 


Re-arranging  1.4  gives 


-  ip  =  cy  +  cuy  +  xiip  +  -£~-  uy  . 

ck  G 


Differentiating  I . 10  and  collecting  terms  gives 


ip 

^2  __  22,  244  ___  2 

D  -2Dpo)U  +  p  0)  u  -  3DEc 

2  2  2 
-  3Ec  po)  u 

I  i 

2  „  2  2  3  ~ 

i  i  i 

2  2  2 

2  4" 

=  ip 

5DEc  u  -  3Ec  pas  u 

+  ip 

DEc  u  -  Ec  po)  u 

-  Dc 


2  2  11  2 
(d  -  po)  u  )  y  +  2pa)  uy 


~  (I . 11)  gives 

2 

1  1  '  1,1  '  1  1  no)  1 

-ip  =  2cy  4-  cuy  +  u ip  +  ip  +  .  (uy  +  y) 

ck  G 


du 


2  (1.11)  gives,  after  re-arranging  the  terms 

2 

II  III  TIT  I  I  1  1  1 

-u ip  =  3cy  +  cuy  +3 ip  +  — r  (uy  +  2y  ) 

ck  G 


(1.9) 


(1.10) 


(1.11) 


(1.12) 


(1.13) 


(1.14) 


Substitute  1.13  and  1.14  into  1.12. 


. 


.  C  ;  .  ,  •  '  . . v< 

\  H  j  -  i  'v0  [Gc] 

\\r  -  - 

c 

{'  ‘  1  .  ;  - •  V 
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■DEc' 


O  ry  »  '  '  '  '  1  9 

u  y  +  5uy  +  4y  +  — j—  (2y  +  4uy  +  u  y 

C  k  G 


3  2  3 

+  Ec  poD  u 


UyIV  +  3y"'  +-f  V-  (uy"  +  2y ' ) 

c  k  G 


Dc 


(D  -  pco2u2)  y*  '  +  2pa)2uy 


/ 


irx  2  2v  2  ,  2  OT_  2  2  2 

(D  ~  po)  u  )  +  DEc  -  3Ec  pa)  u 


(1.15) 


Substitute  X.ll  and  I ,15  into  (B)  to  get 


3  u^  +  u  ijj 


p  o  8  e 

3cu  y  +  3cuy  + 


3  oo)  2 

~'~i  u  y 

ck  G 


+  2u‘ 


■DEc' 


0  0  0 


2  IV  ,  ,  ,  ,  , 

u  y  +  5uy  +  4y  +  _  , 

c  k  G 


(u2y ' ’+  4uy ”  +  2y) 


_  3  2  3 

+  Ec  po)  u 


1  11  •  r\n\  *  '  ! 

uy1  +  3y  +  (uy  +  2y  ) 

c  k  G 


Dc 


(D  -  po)2u2)  y  +  2po)2uy 


/ 


.  2  2,2  ,  2  2  2  2 

(D  -  po)  u  )  +  DEc  -  3Ec  pa)  u 


(1.16) 


which  is  free  from  terms  involving  ip  and  its  derivatives. 


To  get  an  equation  for  y  in  terms  of  u,  substitute  1,16  into 


1.8  and  collect  the  terms.  After  dividing  through  by  the 

2  2 

common  factor  (D  -  3po>  u  )  the  equation  becomes, 


' 


era  10J,  fu  .  it..  .1"  i  j;:/;p9  rxe  :tsp  oT 
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,4 
d  y 

,  4 
du 


Ec3u2(D 


75  5 

pm  u  "  +  Ec*") 


,3 

+  ±JL 
3 

du 


,2 

+ 

-r  j 

du 


+ 


du 


2  2  ,  „„  2, 


Ec  u(6D  ~  4pm  u“  +  6Ec  ) 


(D  +  Ec2)  (6Ec3  +  cpm2u2)  + 


5  4  4 

Ecomu  244 

“—^-7—  -  cp  m  u 

k  G 


2  2 


(D  +  Ec  ) 


k  G 


(D  +  Ec2)  (3cpa)2u  +  -  2M£SUtLJil 

k  G  k  G 


2  4  3 
cp  m  u 


+  Y 


D  4  Ec  ‘ 
ck  G 


3  6  4 


2  2  ^2,^  2  4  2,  0  m  u 

~)(3Ec  pm  -  Dpm  +  2p  m  u  )  -  Js!—r~ 

/  ck  G 


0 


When  u  is  replaced  by  (1-(3X) ,  the  equation  of  motion  for 
a  linearly  tapered  wedge  with  shearing  force  and  rotatory 
inertia  becomes 

1  ^E(i_pX)2(D  -  pm2  (l“pX)  2  +  Ec2)  ) 

cl  dX  V  7 

3  s 

L  iLX  ^E(l-pX)  ( 6D  -  4pm2  (1-px)  3  +  6Ec2)  ) 

tJ  dX  k 


+ 


1  djy 

2  2 
cl  dX 


^(D  +  Ec2)  (6Ec2  +  pm2 ( 1-px) 2  )  -  p2m4 (1- 


•px) 


+  Epm2  ( l-px)  2  -  ( 1-|3X)  4 


k  G 


(1.17) 


k  G 


' 

¥  S  £ 


( 

f-  + 
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3Ep20)4  ( 1-[3X)  3 


k  G 


0 


+  y 


c 


^JL±JiLi iS\  f3  Ec  2  pa)2 

V  ck  G  2  \ 


ck  G 


0 


(l-px)4) 


0, 


(I- 18) 


where  X  *  ~r  . 


APPENDIX  I-B 

APPLICATION  OF  GALERKIN ' S  METHOD  TO  THE  EQUATION  OF  MOTION 
FOR  A  WEDGE  WHEN  SHEARING  FORCE  AND  ROTATORY  INERTIA  ARE 

CONSIDERED c 

By  multiplying  Equation  X0X8  by  c  and  collecting  the 
coefficients  of  the  powers  of  a)  the  equation  of  motion  becomes 


Ax(l-px)4  y  o>6 


+  [a2(1-PX)4  y"  +  A3(l-px)3  y'  +  A4(l-px)2  yl  oo 


4 


4  TV  T  »  •  »  o  «  » 

+  A5(l— PX)  yiv  +  A6(l— px)  y  +  A?(l~pX)  y  +  Ag(l-px)  y 


9  V  9  TV  a  8  8  8  8  ~ 

+  Agy  O)  +  A10  (1-pX)  yiv  +  Au  (1-flX)  y  +  A12y  =  0  (2.3.13) 


.  ( J'  4-  (x^i)  S uq^cSa  +  a)  ;.;  ~  - 
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where  primes  denote  differentiation  with  respect  to  X, 


X  =  •£  and 


A. 


3 

k '  G 


n2  ,  Ep . 

k  G 

■  t2 


3Ep' 

vk  G 


+  p‘ 


l 


A, 


A 


8 


A, 


*  (p  + 

v  k  G7 

t2 

-cA  (^3p  4- 

V  k  G 


“4”  f 3Ec2p  -  Dp 
k  G  V 


O  2A 
B 

k  G 


11Q 


EA 

4 

l 


l4 


A 


11 


6EAc 

P 


4Ecp 

t3 


A 


12 


6Ec2A 

t2 


=  D  +  Ec‘ 


When  the  one  term  trial  solution  was  used  in  the  appli¬ 
cation  of  Galerkin's  method  to  Equation  2,3.13  the  equation 


4>  .  dX 

i 


was  formed.  Since 


0 


( 


Yl(X)  =  a1(6X2  -  4X3  +  X4) 

4>i(X)  =  (6X2  -  4X3  +  X4)  , 


1.19) 


Vi  +  a 
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and  the  derivatives  of  y  were  taken  with  respect  to  X 
Equation  I .19  became 


J  A^l-pX)4  ax( 


6X2  -  4X3  +  X4)  a)6 


+ 


A2(l-pX)4  12a1  (1-2X4-X2)  +  A3(l-px)3  a±  ( 12X-12X2+4X3) 


+  A4  ( 1“(3X)  2  a x  ( 6X2 ” 4X3+X4 ) 


CD 


+ 


A5  ( 1=(3X)  4  24ax  +  A6(l~px)3  (-24a1)(l-X)  +  A?(l-px)2  12a1  ( 1-2X+X2; 


+  A  (1-PX)  ax (12X“12X2+4X3)  +  Ag  a± ( 6X2-4X3+X4) 


CD 


A10(1-PX)^  24ax  +  Ai;l(1-PX)  (-24a±)  (1-X) 


+  A12  12a1(l-2X+X  ) 


(6X2~4X3+X4)  dX  =  0 


(1.20) 


From  Equation  1.20  it  can  be  seen  that  the  constant  a^  is 

a  common  factor  which  can  be  cancelled,  leaving  a  cubic 
.  2 

equation  m  co  .  When  the  roots  of  Equation  I„20  were  found, 
only  one  root  turned  out  to  be  greater  than  zero,  and  that 
value  was  taken  to  be  the  square  of  the  fundamental  frequency 


of  vibration. 
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